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Chapter 1 
INTRODUCTION 


-игрозе 

Uncertainty is the father of the art of decision making. 
As trite as this phrase may seem, it is certainly true today 
with the rapidly increasing tempo of business and the intro- 
duction of more and more complex variables. Truly, "the 
business executive is by profession a decision maker." 

Since the turn of the century, particularly during the 
period following World War II, an awareness of the ultimate 
ineffectiveness of purely intuitive conclusions has permeated 
the research projects in the fields of management. The 
dominant purpose of this continuing research has been the 
elucidation of a rational theory for purposive action in 
Situations of risk and uncertainty. Several methods or 
techniques have been developed to offer a basis upon which 
to compare and evaluate the expected value or worth of the 
separate sets of outcomes under consideration. Using these 
methods, the optimal decision would be the one with the 
greatest expected value or worth. 

To exemplify the methods cited above, consider this 
Simple case. The decision maker has two separate possibilities 
and he knows the probability associated with each possibility. 
The sum of these separate possibilities, weighted by the 
probability that each will occur, is called the expected 


Value of the probability distribution. If the decision maker 


lneferences listed at end. 





has an expected nonetary return for each of the separate 
possibilities, he can determine the "welghted average" of 
each monetary result. Thus he can choose the act with the 
highest weighted average monetary result. This is called 
"maximizing expected monetary value." 

One vital element has been omitted, or tacitly ignored, 
in most of these techniques, namely, the value of the outcome 
for the ທ າ > person based upon his particular funds, 
goals and risk preference. The utility theory has recently 
become accredited aS a means of incorporating the element of 
risk into decision making. 

Various definitions have been ascribed to this theory 
by the different theorists. A utility is what the researcher 
takes to be the value of an outcome for a person. Utility is 
the power to satisfy human wants. It is a measurement of the 
degree to which satisfaction is obtained or it is a measure- 
ment of the degree of pleasure or desire for one object as 
compared to another. Lastly, it is the actual value of a 
dollar to the individual when evaluated in terms of his 
financial status and risk preferences. We shall consider 
utility as synonymous with preferability; thus the degree of 
utility will be taken as the degree of preferability. If a 
person prefers an orange, for example, to an apple, this 
will indicate that he has a higher utility for an orange 
than for an apple. 

The loss or the gain of a dollar has a far different 


Significance to a business man of limited capital as compared 
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to a large chain store concern with unlimited resources. In 
addition, there is a different significance to a man who buys 
only "blue-chip" stock as compared to the man who deals 
primarily with speculative stock. The straight maximization 
of "expected monetary value" theory would not be a satis- 
factory criterion for a decision maker in either of the above 
examples. Now if we have a means of assigning a single 
utility to each outcome, and a Single probability to each 
action-outcome combination, then a decision maker's choice 
of an action can be predicted by maximizing expected utility. 
Thus he has included the element of risk into the evaluation. 

Inherent in the utility theory approach has been the 
difficulty in measuring or identifying the utility value. 

The methods, which have been documented, either involve 
complicated axioms or lengthy laboratory experimentation. A 
completely general method of measuring utility, which can be 
applied in all situations, has yet to be developed. 

The purpose of the work reported in this thesis was to 
investigate ways to make utility calculations by digital 
computer possible. In order to do so it was necessary to 
devise a basic mathematical function which would work for 
all, or at least a majority, of the empirical utility curves. 
This accomplishment could prove of two-fold usefulness. 
First, it would make available to a decision maker an easily 
administered procedure to m n data for individual 
utility curves. Secondly, it would provide a mechanized way 
to compute the utility-dollar values, to determine the 


utility function and to plot the resultant curve, 





ມ 


It 1s hoped that this thesis can lend support to utility 
theory as an aid to decision making and as a rule for guiding 
judicious behavior in the face of uncertainty. If so, possibly 
it will find general acceptance as an economic evaluation 


technique in the future. 


Contemporary Yardsticks 


Payout time is the time required for the cumulative net 
earnings from an investment to equal the investment. Probably 
this technique is the criterion by which the majority of all 
properties are evaluated even though it is inadequate in 
many respects. The primary weakness of this procedure is 
that it considers only the time of recovery and does not take 
into consideration the expected total earnings of the invest- 
ment which depend upon the rate of earnings and the life of 
the investment. Payout time has the advantage of being simple 
to understand and, therefore, is more widely known by decision 
makers. 

Profit-to-investment is the ratio of the total return 
from the investment to the amount of the investment. “his 
procedure considers the total return and thus projects a 
ready balance sheet estimate for the project, but does not 
reflect the rate of return. 

The net present value refers to all future net revenues 
before taxes, discounted back to the date of evaluation. 

This procedure incorporates another valuable dimension into 
the analysis, namely, the present value of annuity earnings 


in the future. Іп the long form approach, the net revenue 





5 
for each year is discounted back to the present, Although 


this step-by-step procedure is more time consuming, it offers 
the advantage of giving a detailed picture of the future 
earnings. The short-cut approach utilizes the basic expression 
for exponential (constant percentage) decline. 

A derivative of the net present value method is the 
profitability index which is the present value of earnings 
divided by the investment. The earnings are discounted over 
the period of the annuities and calculated on the basis of 
the present worth values of such annuities. Some disadvantages 
of this method, which probably have restricted its popularity, 
are: (1) the method requires that an earning rate be 
established before calculations are made, and (2) since the 
index is merely an abstract number, it is not as readily 
appreciated by management and stockholders. 

The rate of return is essentially the interest rate for 
which the "present worth" of the future income is equal to 


the investment. This rate is determined from the formula: 


P = C чту 


1(1 + 


where: i = interest rate for a given period 
n = number of interest periods 
P = principal sum invested at the present effective 
date of the valuation 
C = one of the Series of equal payments made at the 


end of each interest period. 
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The present value techniques have gained recognition as 
better criteria than the payout or profit-to-investment 
because they include the effect of time on the value of 
money. Unfortunately, they do:not include the element of 
risk in the comparison and evaluation processes. Whether or 
not the decision maker consciously includes it in his evalua- 
tion, this element is ever present. Any technique which will 
introduce risk into the evaluation process will be a step | 


toward providing a better standard for the decision maker. 


The Problem 

iuc business terminology includes qualitative and 
quantitative outcomes resulting from decisions made under 
certainty, made under risk, and made under uncertainty. 
Although definitions of these terms are generally consistent, 
the author believes that more clarity will result by specify- 
ing the definitions as used in this paper. 

if each action is known to lead invariably to a specific 
outcome, we have certainty. If each action leads to one of 
a set of possible outcomes and each outcome occurs with a 
known probability, we have risk. If each action that may be 
chosen is identified with a distribution of potential out- 
comes rather than a unique outcome, we have uncertainty. To] 
the outcome is one which, following the choice of a course 
of action, is either obtained or not obtained, we have a 
qualitative outcome. 11 the outcome is one which is or is 


not obtained in various degrees but the extent to which such 





7 
an outcome is obtained is potentially measurable, we have a 
quantitative outcome. 

The simplest type of problem situation is one which would 
involve only two possible outcomes, both qualitatively defined, 
and two courses of action. It is necessary to know the rela- 
tive values of these outcomes to the decision maker in order 
for him to determine which of these two courses of action 
would be better. Recognition of the necessity for determin- 
ing the relative importance of outcomes to decision makers 
has given rise to the value or utility theory. 

As a generalization, it may be said that insofar as it 
can be measured quantitatively, the over-all objective of a 
business adventure is to maximize the profit on the investment, 
consistent with maintaining a sound financial position.“ . One 
of the most prominent "maximization" analysis of rationality 
1s the hypothesis that the decision maker should maximize 
expected utility or value with respect to his beliefs con- 
cerning the facts of the situation. To perform this maxi- 
mization, he needs to have a subjective probability function 
which measures his degree of belief and a utility function 
which will measure the relative value to him of the various 
possible outcomes of his actions or decisions. 

Thus the problem is to derive the subjective probability 
function that will measure these relative values for the 
operator so that he can make the optimal decision which, 
consistent with our definitions, will maximize the mathema- 


tical expectation of value or utility. 





The Hypothetical Model 


Uncertainty in management decisions generally can be 
contributed to at least two major sources: imperfect fore- 
Sight and human inability to solve complex problems contain- . 
ing a host of variadles even when an optimum is definable. 
Under our definition of uncertainty, each action which may 
be chosen is identified with a distribution of potential 
outcomes rather than with a unique one. To help define the 
boundaries of the potential outcomes, many companies have 
‘made use of models. This has permitted the substitution of 
empirical methods for the former purely intuitive or ratto on 
grounds for selecting a criterion. 

A model, by the simplest definition, is a representation 
of reality which attempts to explain the behavior of some 
aspect of it. Because & model is an explicit representation 
Of reality it is always less complex than the reality itself, 
but it has to be sufficiently complete to approximate those 
aspects of reality which are being investigated. The model 
is constructed so that the parameters represent as many input 
variables as possible. By varying one or more parameters at 
a time while keeping the remainder constant, the limits to 
the distribution of potential outcomes can be defined. In 
the past, time and human error have been the major deterrents 
to practical use of any necessarily complex model. With tne 
advent of the digital computer, management now has a convenient 
means of obtaining timely and accurate results even with the 


most complex models. 





en proposed that the decision maker assign a 
subjective, or "utility," value to potential dollar conse- 
quences. To accomplish this, the operator must first con- 
struct a utility function or curve for himself. Utilizing 
this function or curve, he can convert dollar consequences 
of act-events into units of measurements, called utiles, 
which reflect what the dollar really means to him. He then 
assigns probabilities of occurrence to each event. The 
weighted average consequence of each act is found by weighting 
the utiles by the probability. He is now able to select the 
act which will give the largest weighted average consequence, 

We shall use the Grayson thesis to a large extent in 
setting up his model. The model will consist of a short 
questionnaire and an associated computer program, The 
questionnaire can easily be administered to management 
personnel or even to prospective employees. The probabilities 
(answers) from the questionnaire will be converted to utiles 
by the computer program. Then the computer will attempt to 
derive an equation to satisfy the dollar-utile data points. 
Using the derived equation, the computer will plot the 
utility function for the participant. From this point, the 
uses of the plotted utility function will be limited only 


by the practical use desired by the person or firm. 


Plan of the Thesis 
The general outline for the thesis will be the develop- 
ment and detailed description of the model which has been 


introduced in this chapter. 
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The evolution of the utility theory will be traced from 
the early Bernoulli "diminishing utility curve" to the current 
period in Chapter II. A resume of the different theories of 
measurability of the utility value and references to the 
major contributions in this field will also be found in this 
Chapter. 

In Chapter III we develop the mathematical equation 
which will be used in processing the data used in the research. 
The various methods of regression analysis which form the 
basis for the computer programs will be briefly explained. 

The development of a questionnaire, which is our sug- 
gested contribution to the problem of measurability of the 
utility value, will be found in Chapter IV. 

A method of simultaneously incorporating the data from 
Part II with those from Part I is developed in Chapter V. 

The practical usefulness of the questionnaire was 
tested by 200 volunteer students from the School of Business, 
University of Kansas. The results of this test form the 
basis of Chapter VI. 





ona ter II 


THE DEVELOPMENT OF THe UTILITY neon 


Historical Perspective 

The concept of Utility Theory is not a recent топос сое 
Kauder credits Kraus? with discoverlng some early discussion 
of the law of diminishing utility in Aristotle's works.’ In 
Ш 20, the mathematician, Daniel pernonll1 сеа сшита, 
with the games of chance, isolated the theory of marginal 
value as applied specifically to the value of money.? He 
theorized that the value of the last item received of a com- 
modity is the marginal UCM A His "diminishing utility 
curve" equation: 


us =n шо C aa le 


is equivalent to the basic equation for the research of this 
paper. 

Adam Smith made the distinction between value in use 
and value in exchange in 1776.7? Jeremy Bentham, through his 
€ culus of values, introduced the concept of utility as a 
numerical magnitude. He also utilized the principle in his 
Suggestion for the measurement of quantities of pleasure and 
pain, 1° and for the creation of the following corollaries: 
(1) gambling is utility-decreasing and, (2) insurance is 
utility-increasing.+? 

1t was left to Heinrick Gossen to publish the fundamental 
principle of marginal utility theory. He concluded that: "a 


person maximizes his utility when he distributes his available 
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money among the various goods so that ne obtains the seme 
amount of satisfaction from the last unit of money spent 
upon each commodity, "+ 

Jevons was probably the first person to introduce utility 
as a positive and integral part of economic philosophy. In 
the 1860's he attempted to reconstruct political economy with 
his speculative calculus of pleasure and pain. “This concept 
later was termed by economic theorists as the preference 
approach to utility. He said "value (was) to be established 
on the basis of labor and the problem of rent, wages, inter- 
est, etc., (and was) to be solved as mathematical functions. "+3 
La Nauze stated that Jevons' works mark the "beginning of 
eeonomics, as distinct from political economy. (His) Wei lies 
theory made economics a 'serious' subject, a technical subject. "1" 

Walras,1? Marshalll® and Rdgeworthl7 concluded that 
utility was a quantity which could be measurable if the 
Operator had sufficient facts. Pareto definitely abandoned 
this approach; he returned to the Jevons concept of a scale 
of ences. 18 

The first careful examination of the measurability of 
the utility function and its relevance to demand theory was 
made by Fisher in 1992. He resolved the measurability problem 
pon the case in wnien tne marginal utilities Of the various 
quantities are independent of one another, 17 

in summation, the evolution and development of the 
utility theory has had periods of virtually universal appeal 


followed by lengthy periods of dormancy . 22 Following 





ша 

Bernoulli's mathematical derivation of the "diminishing utility 
curve" in 1730, little was documented, if any work was actually 
done, until Gossen's publication sparked the pioneering efforts 
of Jevons, Marshall and Walras. Fisher became the first 
American to seriously advocate and demonstrate the applications 
of the theory. Even with the advent of modern methods of publi- 
cation and distribution, another dormant period existed until 
1926. 
Contemporary Perspective 

Bagnar Frisch could well be considered the father of the 
modern surge of research into the application of the theory. 
He introduced his axiomatic method in his 1926 article under 


el This method has 


the heading of “axioms of the Second kind." 
proved itself of great value particularly in geometry and 
probability. While many developments have taken place in 


22,23,24,25 undoubtedly the 


axiomatizing economic theory, 
greatest impetus of recent years has come from the work of 
von Neumann and Morgenstern^? and Магвоһак.27 

Although the comparison of adjacent preferences appears 
во nave its first explicit recognition by Pareto, 24 the modern 
introduction of this concept appears in Frisch's article. Ти 
1934, after the appearance of the Hicks-Allen contribution“? 
the question of adjacent preferences was revived by Lange29 
followed by Samuelson, 24 Weldon, 32 Suppes and Winet, 33 and 
34 


von Neumann and Morgenstern. 


presented Separate articles on 


e 


In 1954, several authors 


an entirely new approach. Hausner32 and Thral129 are credíted 





14 


with establishing the existence of a lexicographic utility 
function from their investiravion of the non-Archimedan be- 
havior under uncertainty.  Debreu made an analogous investiga- 
tion in terms of the classical commodity space, 27 п ее 
Georgescu-Roegens* work added significance to the concept . 2° 
Chipman defined utility, in its most general terms, as "a 
lexicographic ordering, represented by a finite or infinite 
dimensional vector with real components, unique only up to an 
isotone (order-preserving) homogeneous transformation; and 
these vectors (or ‘lexical numbers') are ordered lexicographically 
like decimal numbers or words in a dictionary. Less generally 
stili the components of the lexical numbers may be additive 

(in the Bernoulli case) or multiplicative (in the Frisch case); 
at any rate, unique up to a homogeneous transformation which 


preserves group operations. "7? 


Measurability 

Because there have been so many thorough and learned 
dissertations in the last decade on the various ideologies 
and definitions appertaining to utility theory, it would be 
impractical and useless to P to condense tne major 
contributions into this thesis. Accordingly, detailed 
explanations of the various theories will be omitted; only 
references to selected measurement theories will be included. 

The concept of a measurable utility, i.e., of a real- 
valued function appropriately linear with respect to proba- 
bility distribution, measuring an individual's preference 


ratings, originated with Bernoulli. A completely rigorous 





12 
formulation and treatment of the existence of sucn a utility, 
on the basis of a well-defined set of axioms or postulates, 
however, was unknown until Fisher’s work appeared in 1892. 


This was followed by von Neumann and Morgenstern in their 


! 


0 
Theory of Games. Marschak developed a set of axioms for 


11 


the case of a finite number of sure prospects. Ruben 


extended the Marschak system to tne case of an infinite 
number of such prospects .*2 For a description of von Neumann 


and Morgensternts axioms and three more recent methodologies 


43 


with their associated axioms, see Ackoff, 


Hick and Allen introduced the ordinate theory as opposed 


ny 


to the then existing cardinal ideas, Furtner studies deriv- 


ing the existence of an ordinate utility function from postu- 
lates about preferences have been made by мола” 


46 


and by 


For lengthy discussions on cardinal verses ordinal 


48 


Debreu. 


utility see Láttle, +7 von Neumann and Morgenstern and 


49 


Robertson. The major difference between Hicks-Allen 
utility and classic marginal utility concerns the assumption 
of the orderability of marginal uvcility. 

Armstrong introduced marginal preference? wnile von 
Neumann and Morgenstern spesis of risk index or utility 
index?! as a measure. Davidson and Marschak, > Luce?3 and 


54 


Papandreau studied stochastic choice behavior as a means 
of measurement. 

A logical approach to a discussion of utility is shown 
by Kennedy?” who says, "Thus we shall hold that utility is a 
quantity, i.e., has magnitude; chat то топ вте» ауа 


that two utilities when added together do not yield another 





16 
utility...we shall not attempt to decide whether utility 
should properly be regarded as pleasure, as happiness, as 
Satisfaction, as intensity ог desire GC aS ш 

Fisher postulates that by the use of his axioms) "осуу 
as defined...does not involve the economist in controversy as 
to the laws of the subjective states of pleasure and pain, 
the influence of their anticipation as connected with their 
probabilities, the vexed questions whether they differ in 
quality as well as in intensity and duration, whether duty 


6 


can or cannot exist as a motive independently of pleasure,.." 


Modern Applications of the Utility Theor: 





An exhaustive resume of the various applications to wnich 
utility theory has found acceptance today would be prohibitive 
in this thesis. <A few examples will be cited below to portray 
the magnitude of the spectrum of applications which reflect 
the versatility of the theory. 

Since World War II there has been a tremendous increase 
of interest in the use of scientific methods in solving 
management problems, under such labels as “operations research", 
"nanagement science", etc. Churchman, 27 Ackoff,29 Miller and 
Starr, 22 Thomas and Deemer, °° to mention only a few, have 
included utility theory as a possible scientific method for 
use in decision making. 

Grayson attempted to simplify ee decisions bY oil 
companies by use of utility theory, Ol Kaufman has extended 
Grayson's work in this 21816. 92 

Koopmans asserts that "it is shown that simple postulates 


about the utility function of a consumption program for an 





i 
infinite future logically imply impatience at least for 
certain broad classes of programs, «92 

Tinbergen attempted, with the help of measured economic 
concepts, to indicate the rate of savings, as a function of 
time, which maximized utility over tine. 7 

Samuelson used value assumptions to help him move 
Straight from the utility of the individual to the welfare 
of society. 55 Plguo and Little made use of inter-personal 
utility comparisons to take them from the utility of the 
individual to the utility о? the community and only then 
introduce value assumptions to bring them finally to economic 
welfare, 9° 

Tne UCLA Design Research Program has experimented with 
the application of the utility function in devising a method- 
ology for defining and applying consistent, realistic value 
quantities useful in decision-making at all levels of the 
design process, 57 

Althougn adoption of the utility theory has been extended 
to an ever increasing number of uses, the process is still in 
its embryonic state. Utility theory can incorporate personal 
attitudes toward risk са пс into the evaluation of decisions. 
It should prove to be of substantial assistance in creating 


consistency in decision making and should provide a basis for 


attempts to make Simple business decisions by computers. 





Chapter III 


THE MATHEMATICAL EQUATION 


Background 

As mentioned earlier in Chapter II, Danicl Bernoulli 
was one of the first to present the general idea of intro- 
ducing subjective values of dollars into oxpectation calcu- 
lations. He rationalized that each equal increment of gain 
yielded an advantage which was inversely proportional to the 
individual's wealth, 1.е., 


dU = K (1) 


м 


where dU is the increment of utility resulting from an 
increment dx of wealth and K is a constant. From this 
ше шретоп, he projected tne total utility as a logarivomic 


function of wealth in this manner: 


Ux) = K log = (2) 


where C is the amount of wealth necessary for existence. 
integrating the differential expression, we get his 


"diminishing marginal utility" equation: 


1 


ເ = k log (x) + C (3) 


where the constant is determined by the conditions that, 


68 Thus Vice 


when wealth is at the subsistence level C, U= 0. 

can be seen that the utility of money was represented as a 

linear function of the logarithm of the dollar value. 
Research done by Kaufmann ©9 on the ЕО fFURNCLIOR 


curves derived by Grayson/O indicated that & logarithmic 
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function was not the best Tit for all cases; In toc e 
stated that a cubic or a quadratic function (ір Soremo R 
personal curves better than а logarithmic curve, while ІП 
other cases, a quadratic mignt best Mit tne negative Scere 
(quadrant III) and a cubic might best fit the positive 
section (quadrant I) of the curve. 

Intuitively it would appear that the logarithmic curve 
would more closely match the curve in Figure (1) than any of 
tne other typical named curves. Without discounting the 
contribution of Kaufmann, we proceeded on the assumption 
that in a majority of cases, the logarithmic equation by > 
Bernoulli would give the better fit. 

Accordingly, the following equation was used to commence 
the rescarch: 


ој) = а log ((х()) > ӘЛЕ) (4) 


The constant which appeared in ວເ" equation was in- 
corporated into the basic ее value by the addition 
of the b in the denominator. ‘This has the effect of forcing 
the curve to pass through the point (0,0). ) (With xo ED NEUE 
logarithm of l does equal 0.) There is one requirement; b 
must always be greater than the largest negative value of 
Х(4) in order to eliminate the possibility of having a 
logarithm of a negative number. Subject to this requirement 
we have successfully incorporated the Bernoulli constant 
Within our logarithm expression. 


Having established equation (4) as the research equa- 


tion, we now can look at the methodology for generating the 
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data points, The expected value of any distribution of 


events is: 


n 
E, = ) P. V (55 
мен 


where P, is the probability that the event takes place and Vy 
is the value of that event. In our case where the value is 


measured by utility, the expected utility value is: 


n 
E, = ) PL Uy (6) 
K=1 


where Uy, is the utility of the particular event. 

In an investigation of utility, a useful concept is that 
of indifference. A person is said to be indifferent to 4 
proposal if he does not care whether he accepts it or not. 
me he is indifferent it is axiomatic that the expected value 
of the proposal is zero for him, 

In the case where the proposal has two possible out- 


comes the expected utility is given by: 
Ey = Pj U) + Po Us : (7) 


If a subject is indifferent to the proposal then В) = 0 and 
P4 and Po are said to be indifference probabilities. 
Indifference preferences have been a recognized theory 
of measurement within the utility theory. Grayson utilized 
the indifference probabilities to determine his utility-in- 
money functions. We have used it in both Parts I and II of 


the questionnaire where the expected value of the -hypothetical 
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wager was set equal to zero. Thus: 


Uwinning Pwinning Е Ulosing m winning)? = 0 (8) 


and P ís the indifference probability. 


winning 


General Mathematical Model 


The equation of the.curve is: 
а = anal 95253] (9) 


The objective is to determine a curve such that the sum of 
the squares of the distance of the data points is the least, 


where e = data points. So let: 

ej = uj -a ln [(x, + b)/b] (10) 
Now squaring both a 
e; = u, -aın ((x, + b)/b) F (11) 


Y 
And squaring over all points: 


n 
i % (e) (12) 
jel 


We take the partial derivatives of E with respect to a and b 
and set them equal to zero in order to minimize the total 


error. The steps in sequential order are listed below. 


, 


n 


— у (e) 9-1 20 (13а) 


j=1 


Hs 
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n 
IE _ Qo aan 
к = 2 > (е) = 0 (139) 
j=l 
294 == [(х, + b)/b] (14a) 
0 
w^ = a [1/(x, + b) - 1/0] (14b) 


Substituting equations (11) and (1^) into equation (13) and 
dividing by -2, we get: 

n n 

) w/Ggeb)-1/e]-a ) InU(xj+b)/0] [1/(x,+0)-1/0] = 0 (15а) 
4-1 j=1 


n 


n 
) u, [1n ((x, + b)/b)] - a ) [1n (x, T b/b) ]? = 0 (15b) 
j=1 B 


To simplify the equations and to reduce repetitive, calcula- 
tions by the computer, the following terms were collected 


and identified as follows: 
Sam, = In [(x, + b)/b] 
H, = [1/(х, * b) - 1/9] 


Now rewriting the partial differentials after introduction of 


the collective terms, gives: 
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n n 
ЗЕ = > uy Hy ~a > Sam, H; = 0 , (16a) 
TES a 
n n | 
IE _ _ + ,. | 
SCR > u; Sam, a > (Sam, ) = 0 (16b) 
gm J=1 


In order conveniently to prepare subprograms for the use of 
the mainline computer program, the following function 


designations were used: 


j=1 
n 

Y ща ^ u; H; 
3-1 

n 
j=1 


Rewriting the partial derivatives after the introduction of 


the function designations, gives: 


ЗЕ - т - аВ - 0. (17a) 
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SË = а-ай-0 = (175) 

and 
a = T/R (18) 

ог 
a= Q/W (19) 


so the first derivative of the remaining unknown constant, 
BD, 13: 
FOFB = (T W) - (Q R). (20) 


Тһе Newton-Haphson Method 


When the derivative of f(x) is a simple expression and 
easily found, the real roots of f(x) = O can be computed by 
the Newton-Raphson Method. It is a powerful method of solv- 
ing an implicit equation. This method was used in the computer 
program which has been identified in Зи este as the 


"Original Logarithm Program," 


Y 
The Twenty Question Method 


In this method of iteration, two numbers xy and x, are 
found so that the root of the equation lies between them. 
Since the root lies between the two numbers, the graph of y 
= f (x) must cross the x-axís between x - X4 and x = X», and 
E and yo must have opposite signs. Having determined X4 and 


Xo, we find a mean value for the two points by: 
\ 
хы = (ху + х2)/2. 


We substitute x, for x in the equation y = f(x) and calculate 
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a new value for y. For example, let us assume y, was negative 
and. yo was positive. If Це is negative the value of Zr is 
Substituted for x] and the procedure for finding a new да 15 
repeated; if y, is positive tne value of x, is substituted 
for x and the procedure for finding Xm iS repeated, etc. 
Equation (20) of our mathematical model is the general 
equation with the constant B as the only unknown. To apply 
tnis method, FOFB plays the part of y and B plays the part of 
x. The first Steps are to determine two values for B such 
that FOFB of one (BG) is larger than zero and FOFB of the 
other (BL) is less than zero. Having established these two 
limiting points, we used 20 iterations to calculate the cor- 


rect value of B. If the absolute value of FOFB was less than, 


.00001, the iterations were terminated. 


Snolse Method 


A computer program designed to solve non-linear simul- 


71 


taneous equations was developed by McLean. This subprogram 


was used with the mathematical model for this thesis. 


S 
e 





Chapter Ty 
THE QUESTIONNAIRE 


Very little has been documented on the exact methods 
proposed and used in determining the utility function of 
individuals by other authors or research personnel. Grayson 
used a table of indifference probabilities on hypothetical 
2 


ari l Ling propositions. Mosteller and Nagee made extensive 
laboratory experiments with selected groups of students and 
National Guardsmen. The experiment consisted of a gaming 
technique with small amounts of money (from 1 to 100 cents) , 
designed to determine a numerical value of utility for the 
different sums of money. After obtaining for each individual 
a utility measure for various amounts of money, they attempted 
to predict how each individual would choose among a set of 
uncertain prospects, where the entities were amounts of money 
with associated probabilities. 2 Although several measure- 
ment techniques were mentioned in the previous chapter, none 
of these theoretical techniques were transeribed into simple 
usable experiments or questionnaires. 

The questionnaire used for the computer program in this 
thesis is included as Appendix (1). It was designed as a 
series of questions about investments. To simplify the 
Situation, however, two slightly different games of chance 
were substituted for investments. | 

For Part I, the questions concerned a ກາ game 
played with a deck of 100 black or red cards. The results 
of this game depended upon one and only one card being drawn 


by the hypothetical opponent. Should the opponent draw a 
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black card, the participant won; should ne draw a0 rodgedrgP 
tne participant lost. The answers to the series of entries 
represented the odds or probabilities requested by the par- 
ticipant. The answers reflected the exact conditions the 
participant would require in order to make a monetary invest- 
ment in the amount of the gaming рер от the entry 

For Part II, the game consisted of a series of bets on 
one toss of a fair coin. In this part, the amount of loss 
was fixed and the participant was asked to write the amount 
of win he would require in order to play the game or make an 
equal monetary investment with the probability O1 50/90 01 
winning. Р 

Although we shall defer the detailed explanation of 
fare point until later in the thesis, it should be noted 
that the calculations including only Part I data points are 
much easier than the calculations including Part II data 
points or a combination of Part I and Part II points. In 
fact, the calculations of Part II depend upon Part I results 
for reference. However, initially it was not known which 
method would create for the participant the more meaningful 
game and thus give the more ‘accurate personal utility curve 
` so Part II points were included not only for cross reference 
purposes but also for analysis purposes. 

The basic development of the theories underlying the 
questionnaire can best be explained with reference to Figure 
(1). This figure shows the order of determining the values 


of the points from Part I of the questionnaire. It is to be 
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noted that point (0,0) vas assigned as was point (410,41). 
It snould also be noted that the next points are determined 
in this order, -10, +20, 44.“ освен 

The points (0,0) and (+10,+1) were arbitrarily set. As 
the utility of zero dollars logically would be zero, the 
criteria for the establishment of the zeroth point need not 
be explained further. The dollar value for the first point, 
to represent а, utility of +1, could be any arbitrary amount. 
However, several factors must be considered in picking this 
point. The number of entries on the questionnaire should be 
kept to a minimum; a long questionnaire immediately creates 
an undesirable effect upon the participant. Yet it is es- 
sential to obtain a reasonable number of points in order 
properly to represent the curve. It is recommended that a 
number between 12 and 20 points be used. Lastly, to make 
the questionnaire applicable to persons of varying financial 
Status, the largest monetary entry must be restricted in 
order to keep its value within the rational understanding of 
the С. Thus if a small number like $1 is picked, 
the questionnaire is either too long or the range of dollar 
values is ee for practical use, If a large 
number like $50 is picked, the larger monetary limitation 
can easily be exceeded. | 

Having established the first two points, the other 
points are determined from probabilities and dollar values 
associated with a previously determined dollar-utile point. 


The equation used in determining the other points was the 
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indifference equation: 


Uwin Pwin 7 Úlose (F = Prin’ =? (8) 


Inherent in this equation is the restriction that thg 
utility of one of the points must be known. Thus each 
questionnaire entry must contain a dollar value for which 
the utility value is Know, The logical progressive steps 
for the entries would be; 

(1) compare winning $10 (the initially established 

point) with losing plo: 

(2) compare winning 120 with losing $10. 


(3) compare winning $20 with losing $20. 
(4) compare winning $40 or $50 with losing $20, etc. 


This technique was used in Part 1 of the questiommaires 
however, intuitively, it was reasoned that by representing 
the entries in such a progressive step procedure ШЕ would 
encourake the participant to establish the first probability 
and merely increase the probabilities by some predetermined 
increment. This would tend to defeat the purpose of having 
each entry determined solely on the basis of that entry 
alone. The order of the entries were made random to minimize 
this purely mathematical approach. 

By application of Equation (8) through the logical 
progressive steps, all values of points from Part I can be 
calculated starting with the,arbitrarily established point, 
which we have set as (410,41). Hereafter this procedure 


will be called the step-by-step procedure. 
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Referring again to the equation and Figure (1), it 
Should be noted that each ERA step ís based upon: 
(1) the calculated utility value for the previous step, and 
(2) the probability for the corresponding entry on the ques- 
tionnaire. Thus, the validity of the points and the ultimate 
validity of the derived function equation depend upon the 
validity of each questionnaire entry. Part II was designed 
to provide a reference for comparing the probabilities and/or 
results of Part I. By allowing the participant to establish 
the monetary value of the win, there was no way ог computing 
utiles step by step as was possible in Part I. Therefore, it 
was necessary to have a predetermined Scale of values for 
reference. Part I furnished the only available scale source 
so the utility for each play on Part II was set equal to the 
utility reflected on Part I. Thus it can be Seen that by 


our construction of the Part II questionnaire, 


UX1(1) = - UN) (21) 
ວ. 


ບ шт (2la) 


where UXl; is the utility value for the Part II point and UN y 

is the value of the negative Part i point.  Hereafter this 

procedure will be called the direct substitution method. 
Figure (2) shows the calculated points of Part II as 

(0) points and calculated points of Part I as (A) points. 

Although Figure (2) reflects only the results of one ques- 


tionnaire picked at random, the majority of the questionnaires 
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followed tre example in that Part II gave a curve Walch ae 
larger utility values for corresponding dollar values than 
did Part I. 

AT his point, sovoral alternatives Mere Терс с 
соата use only Part I; c could usc only Part ІЛ” ос ме есіне 
use some method of data smoothing to combine ‘the results of 
Part I and Part II. Although conclusions were made on the 
comparison of the results of altcrnatives (1) and (2), it 
was deemed advantageous to utilize alternative (3) as the 
most logical method of analysis for reasons stated below. 

It was noted earlier in this chapter that two aistinet 
propositions were includes in the questionnaire in order to 
provide a broader base for calculations. In order to retain 
the usefulness of having the two different typos of questions, 
care must be taken to equalize, as much as possible, the 
weight given to each type in the final analysis. In the 
absence of full equality, at least the operator should know 
tne weight given to each. 

The equation for calculation of Part I data points gives 
eompdlete weight to the Part T technique. The direct corrola= 
tion of Part II points to the negative Part I points probably 
incorporates 50% weight from Part I and 50% weight from Part 
II. If the data points were determined for all of Part I 
and then the data points for Part II were calculated from 
the results of Part I, any data smoothing process for these 
points would result in a higher weighting factor for Part I 


than Part II. To overcome this weighting factor completely 
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is impossible since it is necessary to use the values of Part 
I as a reference scale in determining Part II. Therefore any 
method which would introduce individual points of Part II, 
during a Part I type calculation would reduce the weighting 
factor of Part I. 

A possible method of solution would be one which would 
calculate the dollar-utile points for Part I progressively 
from zero to the dollar point just immediately larger than 
the smallest dollar point for Part II. At this time, by 
means of a data smoothing technique, include the Part II 
point and recalculate all points on Part I to correct for 
thls one point from Part II. Then the normal step-by-step 
method of calculation of Part I points would be resumed 
until the next larger Part II point was past. Recalculating 
all of the previous points before proceding would insure a 
more equalized weight from Part II. 

Again several alternatives were available for projecting 
the points of the Part II curve into the final curve. An 
equation could have been determined for the Part I points; 
likewise, an equation could have been determined for the 
Part II points. Then a mean could have been found for all 
the points from these two equations. Several reasons rendered 
this solution impractical. If the equation included only the 
positive points of Part II, insufficient points меге avail- 
able to make this procedure practical. If the positive 
points of Part II and the negative points of Part I were 


used to determine the curve, the added weight factor for 
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Part I is still present. With this latter method, modifica- 
tion of the combined curve as well as the future Part II 
points after each separate point of Part II was incorporated 
into the combined curve would have been virtually unrealistic 
if done on a step-by-step basis. Several methods for calcu- 


lating the final curve are explained below. 


Part I Alternative 

The simplest method would be to analyse just the Part I 
points. The dollar-utile points could quickly be determined 
by using the step-by-step procedure. The final utility curve 
could be determined by use of the various methods of regres- 
sion analysis discussed in Chapter III, This alternative would 
be applicable for the participant who places more reliance on 
his responsiveness to the type of entries found on Part I or 


uses this type of reasoning in his daily business decisions, 


Part II Alternative 

As the name implies, this alternative utilizes the data 
from Part II. Part I points can be calculated by the step- 
by-step procedure. Then Part II points are calculated by the 
direct substitution procedure. To be useful, a curve must 
have both negative and positive points; therefore, the only 
effective way of developing this alternative would be to use 
the Part I negative points and the Part II positive points. 
This alternative is applicable for the participant who places 
more reliance on his responsiveness to the toss of a coin 
technique. It would also be applicable for the participant 


who normally makes business decisions of this type. 
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Part III Alternative 

With this alternative, all of the points of Part II are 
incorporated with the Part I points. For the operator who 
does not have sufficient information to judge the participant's 
psychological reaction to the two different types of entries, 
this alternative would be most effective. Likewise in the 
case of a participant who places no added reliance on either 
method, this alternative would be most effective. 

We tried two methods of data smoothing for this alterna- 
tive. Both methods can be explained more easily by referring 
to Figure (3). 

Tne rirst method is based upon the normal step-by-step 
calculation method for points of Part I with the periodic 
incorporation of the points for Part II. The normal step-by- 
step procedure for calculating Part I is started. After each 
point is calculated, the dollar value of that point is 
compared to tne dollar value of the first Part II point. If 
the Part II point is larger, continue the normal procedure. 
When a Part I point is found to be greater than the first 
Part II point the normal procedure is terminated temporarily 
to incorporate the Part II point. Considering the incorpora- 
Zion of the first Part Il point (UX14) on Figure (3), the 
function for Part II was supposed to be linear from point 
Ко through point UXl4 and projected to X4'. Ihe mean value 
of X4 was calculated using the projected value ໃຊ and the 
normally calculated Хз. Because Хз" depends upon СА» this 


latter value was recalculated. Now we are ready to continue 
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with the normal step-by-step procedure until the next Part II 
point is reached when the same recalculation procedure is 
made. 

All previously calculated points should be recalculated 
whenever a Part II point is incorporated. Herein lies the 
difficulty of this method. To recalculate each previous 
point by this method would be very time consuming. 

All the methods described in this chapter have used a 
Step-by-step procedure for calculations. In the next chapter 
we shall describe in detail a method whereby the utilities 


of all questionnaire data points are computed simultaneously. 
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Figure 3 


STEP-BY-STEP PROCEDURE SHOWN GRAPHICALLY 





Chapter V 
SIMULTANEOUS UTILITY CALCULATION 


Because of the difficulties inherent with the step-by- 
step methods illustrated in the last chapter, we would like 
to develop a procedure which would make stepwise calcula- 
tion unnecessary. We choose to do this by adopting an 
approach whose basis was that of the general least square 
concept. 

In this method, we write a series of equations which 
will include each of the utility values (Us, UN;, 011 ] ) when- 
ever they appear at each data point. By solving these equa- 
tions simultaneously, each value of the variables is considered 
in the calculations thus there is simultaneous calculation of 
each point in the curve-fitting process. 

The preceding matnematical illustration and the discus- 
Sion of the mathematical model have been based on the premise 
that the participant analyses each questionnaire entry with 
true indifference. Unfortunately, зе premise is invalid. 
For this reason Equation (8) is not meaningful. We must 
incorporate the contribution of this non-indifference 
equation, 


We start with; 
Uwin Pwin 7 Yrose (1 - Prin? 


If we had a true indifference situation, the function would 


be set equal to zero and we would have Equation (8). However, 
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because we do not have true indifference, we set it equal to 
ejs which represents the indifference error. Although the 
value of ej is not known, it is hoped that it will be as 
Close to zero as possible. The mathematical notation j 
denotes the jth equation in a series of equations, so we 
can write: 


U; P; - Us (1 - P4) = E . (8a) 


where ປງ is the utility value of the point and P; is the 
probability of winning. 

It was noted in the discussion of the step-by-step 
method that one of the utility values must be known in order ' 
to use Equation (8). From this it can be seen that the 
Calculation process alternates between the use of a known 
positive value and the use of a known negative value. To 
make the calculations easier to follow, U, has been designated 


as the utility of a positive point and PR, as tne probability 


j 


of winning in the case where we proceed from a known positive 


U, to determine the value of a negative utile. (These PR 


J 
probabilities result from the non-equal money entries on Part 
I.) For the reverse procedure of calculating a positive Us 
from a known negative value, UN, has been designated as the 
utility value of the negative point and P; as the probability 


when calculating an unknown positive value. (These P, proba- 


J 
bilities result from the equal money entries on Part I.) With 
the value for UN being negative, we can now write the equa- 


tion for calculating a positive value as: 





0 


and the equation when calculating a negative value as: 


U, PR, + UN; (1 - PR.) = е, . ( 80 ) 


"E 


In the discussion of the direct substitution method 
for calculating the Part II points, we suggested the 
equations 


prese = - UN (21a) 


where UXL, is the utility value for the jth Part II point 
and UN; is the utility value for the jth negative point from 


Part I. Incorporating the indifference error, we get: 


We have assigned a value to ປງ = 1,0 and called it the 
zeroth variable, so Uz becomes the first variable or unknown. 
In order to establish the equations for the least square 
method, we will group together all of the equations of a 
type, for example, the equations used to determine a negative 
value from a known positive value. If we assume n number of 
data points for Part I and m number of data points for Part 
ІІ, we can write our series of equations beginning with the 


first positive term, U5, as follows: 
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es = U4 Рз + UND (L - РЗ ) (B2 ) 
Cn-1 = Un Р + UN, -1 (1 = Pa ) (B, zl ) 


Con-1 = Uy, PR, + UN, (1 - PRA) (Bon-1) 


Con = UXT + UN] (Boy) 


Congm = UX + Un (B2n+m? 


The above set of equations does not represent the complete 
series for this method. The equation which will couple the 
UX1L's to their neighboring U's will be included following a 
detailed explanation of the above equations. It is recognized 
that we have a sufficient number of equations above to solve 
for the unknowns by setting each equal to zero. However, 
use of е, becomes a necessity when the coupling equations are 
added to this set. 

Following the procedure described in the general 


mathematical model illustration, in the previous chapter, 
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the summation of the squares of the indifference errors: 


LL 


ae 


Es ) ( n 
= 2 
J=1 
where LL equals 2n + m - l. 
We now take the partial derivatives and determine at 


what points these derivatives are not equal to zero. 


LL 

(a) ee =) = from k s 2,n 
j=l 

де 


here — = 0, unless: 


> де 
| = ni Ei 
ог: 
= k = 1 wh ш Р 
j = mn + mn В 
LL 
ОЕ «> е e 
vr 2 ອງ ໆສິ: from k = l,n 
j=l 
OE 


here SUN, = 0, unless: 


0 
j = K where pl = (1 - Р,,1) 


u де 
| ва + К = (2 - PN,) 


ог: 
де 


s 2n 4 k - l where ——À s 1. 
; SUN, 





ЇЇ +3 


де 
E) : В - 
(6) SUXIL 2 e, рапа from k = 1, nxl 
j=l 


де 


here са = О, unless: 


k 
де 
= 2 k =- h = l. 
j n + К l where sol al 


We have a total of LL number of equations in our 


illustration. The first (n-1) equations come from; 
LL 
де 
K 
3-1 


The second (n) equations come from: 


LL 


yom! 
e o 
UN, ^1 
J=1 


| 
O 


The last (m) equations come from; 


LL 


3 saa! е; OF 


j=1 


i 


We now establish an array: 


Ау 111 + А 272 + A, 33° ere ee = ວ 


where: 





LL 


4h 


Ural where k < n 


where n< K< Zn 


k-n+1 


where 2n < К <2п+ш 


Ux. 


(U, is the first unknown because U; = 1.0 is known.) 


To exemplify this procedure, we now look at the array 


positions which receive values from the partial derivative 


Зе 
Б. щ 
jal 


but this is equal to; 
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де де 2 
5 % E SU, елге = 0 





(All other derivatives in the sum are zero.) Note that the 
second equation comes from the derivative with respect to U, 


because U, is the first unknown, 05 the second, and so forth. 
AS 82 = U4 Рз + UND (1 = Рз) 

and en+2 = 03 РЕЗ + UN; (lL - PR4) 

the second equation is: 


2 2 | > 
ШЕН + Рз ) 93 + Pz (1 - P4) UN, + Paz (1 - PR.) UN, = 0. 


3 


0 


50 20551) 
A(2,2) = PR,” + E 


A(2,n) = PR3 (1 - PR.) 


A(2,n+1) = Ра (1 - Pa), etc. 

The next portion of the discussion will involve the 
procedure to be used to incorporate the non-indifference 
equations. It will be noted that there is no way to predict 
wnere the Part II points will lie in relation to the ? 
points; thererore, the same index limits which were used in 
the (b) equations do not apply here. We will use & sub- 
program to search for the relative locations of the Part II 


points. After establishing the location of each Part II 


point as on or between two identified Part I points, we can 





i 
ена э = 
E I => = - 
e 
o е s => 
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Figure 1 
GRAPHIC MODEL FOR THE COUPLING EQUATION 


write an equation which will enable us to add the contribution 
of this Part II point to the immediately associated Part I 
points. In Figure (3), a point from Part II (х1,) lies 
between two points from Part I (ໄຖ, and Хут). We will con- 
sider the function is linear between these three points, 

Then that e which represents the added contribution from the 
Part II point is given by: 
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с, = Uy 25) ບ m - Ud, (22) 


where LJ - X - 1. 

The corollation between the (b) equations and the coupl- 
ing equation can best be understood by reference to Figure 
(3). Xl; will effect the: 

(LJ - 1)th equation, the 

(LJ)th equation and the 

(2n + j - L)th equation, 
where the value of j has no connection with the value of LJ. 
The only knot relationship is that xl is immediately 
greater than X(rj). The contribution to the (LJ - 1)th 


equation is: 


де. 
h nee 
е. (лу) Where (жу = su; 


fe contribution to the (LJ)th equation is: 
2, , 


and the (2N + j =- 1) Ба 15: 


"E 


For illustrative purposes we shall avply the above 
equation to the case where the 2nd X1 (j=2) voint lies 
between X4 end X, (thus ГУ will equal 3) and n equals 5. 
Tne contribution of this e to the (LJ - 1)th equation will 
be: 


n & h 
AX [Up I + Up; (ag) - 924 
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and to the (LJ)th equation: 


ເ f í n Z 
U ບເທ < 


and to the (2N 4 J - 1)th equation: 


-[UNNJ (x) + Uyyg—a (да) - UX1,] 
where (2n + j - 1) = NNJ. 

These separate expressions must be added to the term 
already computed for the three equations. In order to better 
Show the effect of this contribution upon some of the same 
array positions as used by the previous example, we will set 
LJ equal to 3. The following would be added to the (LJ-1)th 


equation in the array positions as indicated: 


0 


I 


Ao 1) 
Alo 2) = Уеа 
Alo 4) = (n/AX)(8/Ax) 


Ao y) = 0 

A(5 3) = -(h/AX), etc. 
and to the (LJ)th equation: 

Al, 7) = 0 

AC, 2? - -(Һһ/лАХ) («/ АХ) 

А(5 3) = (g/AX)* 

Аз д) = 


A(3 11) = -(8/AX) 
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and to the (2n + j-1)st equation: 


(as 

A(3.2) » (AX) (e/AX) 
Alaa 3? = -(g/AX) 
A © 

Ат тъ 1. 


We have only shown the contribution from the non- 
dndi f ference equation., Тһе total value for any от ле arran 
positions will be the array entry previously calculated 
mums the contribution. | 


All other Part II points for which Z, 2 Xl, are in- 


7 
corporated in a similar manner. Once these values are 
determined a library routine is used to solve the simul- 
taneous eguations and compute the utile values for each 
doilar point, 

Thus the calculation of the utility curve for an 
individual can be accomplished by whichever method the 
operator believes will give the more accurate results for 
that individual. Further experimentation may generate suf- 


ficient data which could lead to some qualitative conclusions 


on the use of the several alternatives. 





Chapter VI 
RESULTS OF QUESTIONNAIRE RESEARCH 


The Questionnaire 





Members of the staff and fellow students in the depart- 
ment were asked to complete the questionnaire in the initial 
development period. When the final revised version was adopted, 
about 200 questionnaires were distributed to selected groups 
in the School of Business, University of Kansas. Those 
questionnaires, which were completed and returned by these 
volunteer students, form the basis for this research. 

The purpose of distributing the questionnaire was not 
to obtain statistical data on tne students. We were in- 
terested in testing the usefulness of the queStionnaire as a 
practical means of obtaining data. Our second purpose was 
to investigate whether this data would be suitable for 
computer analysis. 

The scope of this chapter has been limited to the 
questionnaire, itself, and not to the volunteer students. 

We analysed the questionnaires to see if we could identify 
those types of errors which could be corrected in the future. 
Some of the more important types of errors have been documented 
in this chapter. We shall rinish the chapter with a brief 
resume of the success obtained in adapting the data to various 
methods of curve-fitting. | 

In order to have acceptable data for the computer pro- 
grams, a set of minimum standards or limitations were 


established. Those responses or questionnaire entries which 





SE 
do not meet the minimum standards will be identified by type 
of error in the subsequent analysis of the data. 

It should be noted here that the purpose for the ques- 
tionnaire was to create a workable aid for business executives 
to introduce the element of risk into business decisions. We 
shall make frequent reference to this purpose when we discuss 
the types of errors found on the research questionnaires. 

Any questionnaire with less than 4 equal money bets and 
4 non-equal money bets was not analyzed. This number of 
points was considered the absolute minimum required to give 
any realistic results. This type of error was more prevalent 
than any other type. Two participants made the comment that 
based upon their present financial status they could not 
afford to make bets over 210 or $20. Others considered the 
monetary figures unrealistic for a student group. This 
criticism has value, however, it is believed that this 
reasoning reflects, in general, a misunderstanding of the 
value of money. Four bets would involve the maximum loss of 
5100. Although the participants did not have this amount ОҒ 
ready cash to bet, it is believed that any of them would have 
been tempted to borrow this sum in case of loss for the 
expected value of return when the odds were 5 or 10 to l. 
in the cases where the participants are business executives, 
this type of error probably wiil not exist. 

It was surprising to find so large a percentage of 
questionnaires which reflected mere carelessness or a general 


misunderstanding of the Part I type of question. ‘Several 
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cases contained probabilities of winning of 25% in over one- 
half of the entries; these probabilities were noted in the 
lower monetary bets. One of the students indicated in his 
comments that he wondered if we were trying to quiz him on 
his knowledge of the course in probability. This reaction 
was reflected in another case where a probability was shown 
as 28 4/7. 

Another type of error in this general area was the 
isolated entry which was completely unrealistic when con- 
sidered with the other entries. Some of these cases probably 
can be attributed to mere carelessness where the participant 
entered the probability of losing vice winning. For this 
reason alone, it is believed that at least half of the entries 
in this category could have been corrected if the questionnaire 
had been reviewed with the participant when he finished. 

Another problem area resulted when the probabilities for 
“10 values were between 85% and 90%. The resultant curve in 
such a case is unrealistic and would serve no practical purpose 
because it becomes virtually asymptotic to a dollar value of 
very small magnitude. Although this might in fact reflect 
human behavior for a very conservative person, probabilities ' 
of this magnitude appear unrealistic for business Students 
and most likely would not exist in the case of business 
executives. 

A minimum number of participants simply would not make 
an even bet, while one person would make only even bets. In 


each case, the entries on Part II were consistent with the 
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answers reflected on Part I. The present mathematical solu- 
tion for the data from the questionnaire will not work in 
either of these cases. 

Although the proportion of questionnairs which had to 
be rejected was higher ລວກ the response, in 
general, to the questionnaire was satisfactory. The per- 
centage of cases resulting from apparent misunderstanding 
might well reflect on the instructions to the questionnaire. 
In the revision period, the instructions had been redrafted 
several times to incorvorate the constructive criticisms of 
the early subjects. Even with the ເກກ of the 
students, we believe that the instructions contain sufficient, 
easily understood details. Likewise we believe the instruc- 
tions to be of optimum length in that they contain sufficient 
detail and yet are not too long to cause an initial adverse 


reaction from the participant. 


The Curve-Fitting Processes 


The next step in the analysis concerned the curve-fitting 
process on the data from the questionnaires. The use of only 
Part I points proved more successful in the calculation of 
the two constants of our general mathematical equation. This 
had been expected. Hou eni the results even in this area 
Were dependent to a large degree upon the nature of the curve. 
Those questionnaires which contained probabilities close to | 
the even bet situation were the easiest to fit.’ For those 


cases where the later bets were in the 90's, the slope of 
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the curve changed radically and rapidly upward. Unfortunately, 
each of the programs attempted to fit the end points of the 
curve and left the greater deviation in the area where we 
would be concerned the most. By disregarding these latter 
points, the curve-fitting results were much better. 1115 
procedure is not without merit. The sudden increase in 
probability could indicate that the monetary gain or loss 
probably has exceeded the point at which the person has a 
realistic appreciation for the value of that monetary amount. 

The results of the combined curve of Part II points with 
Part I negative points depended generally upon the number of 
points on Part II. As the number of points increased in 
quadrant III, the general tendency was for the amounts to 
increase very rapidly around the 7th to 8th point ($1000 to 
32000). Were the number of points on Part II to equal this 
number, the sharp incline of the curve resulted as it did 
in Part I. Generally the utility obtained for Part II 
points was greater than the utility of the corresponding 
Part I points. This factor apparently made the curve more 
nearly the normal logarithmic curve because it was noted 
that over the same range of dollar values the standard 
deviation was generally less for Part II points than with 
use of Part I points. 

The combination of Part II points with Part I points 
by the step-by-step method created a series of steps in 
both the first and third quadrants. Although the calculated 


curve, using Equation (4), probably was a relatively smooth 
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curve, the standard deviation was much higher than the case 
when only Part I points were used. This step departure from 
the normal curve had its adverse effect on the calculation 

of an equation in each of the programs. The continual change 
of the slope made curve-fitting difficult. 

Of all the programs, the FOFB approach responded better 
than the others.  SNOLSE was dependent, to a large degree, 
upon the accuracy of the initial guesses for the coefficients 
in the case where the curve was erratic. If the deviation 
of the guesses were greater than about 15%, the program tended 
rapidly to increase the value of the coefficients until the 
size of the coefficients exceeded the limitation of 1 x 10415, 
The original logarithmic program was the least responsive of 
all the programs. With calculated data, the program would 
give extremely accurate results if the initial guess for 
BTRY was close to the true B. However, even with calculated 
data when the initial guess for BTRY had a 20 to 25% devia- 
tion from the actual B, the coefficient B was made progressively 
larger as was the case of SNOLSE. 

Our results were indecisive as to the number of coef- 
ficients to be used. We tried to use the equation 


Us = aj ln ( (x; + B)/B) + 22 X4 T a3 (427 


with the SNOLSE program, A majority of cases gave the best 
fit with the use of only the first term containing the 
constants а) апа Б. In a few isolated cases the inclusion 


of the ao term gave better results. The usual result with 
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tne use of the а. terms was that the 84 coefficient would 
increase about two fold before a singular matrix position 

was indicated. With the use of the аз term, two results were 
common, The coefficient 83 became So large or so small that 
its value would exceed the limitation of 1 x 10%) ana give 
us no answer, The second result was a singular matrix. 

In the cases where the curve made rapid ascent in the 
first quadrant, a straight least square program generally 
proved unacceptable. Better results were obtained with the 
use of our basic equation in these cases than with the use of 
a straight least square. The standard deviation was larger 
when we used the basic equation in these cases, but with the 
use of the least square method, we found several instances 
where the standard deviation increased with an increase of 
coefficients. This indicated some error in the application 
of least square to this problem case, 

To insure no misunderstanding on the use of terms, our 
discussion in the previous paragraph was on the general 
least square method of calculating coefficients. Hereafter, 
we Shall analyze the effectiveness of the special method 
which we identified as the least square method of incorporat- 
ing the two sets of data points. As a generality, it can be 
concluded that our least square method was more satisfactory 
than the step-by-step method. Of foremost importance is the 
fact that the calculated data points created a smoother 
curve. This factor made curve-fitting easier: The simul- 


taneous calculation procedure gives equal weight to each 
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value of the variables and thus offers the maximum equality 
to both Part I and Part II points. 

Some difficulties and some unresolved questions were 
encountered with our least square method. In about one- 
fourth of the cases, the method was unstable. The instabilit? 
was indicated by a singular matrix as an answer or a change 
of sign for all of the utile points. It is believed that the 
cause for this difficulty arises from the inability of the 
general least square approach to respond satisfactorily to 
erratic data points. It was noted that in some cases the 
calculated utile values for corresponding dollar points were 
from two to three times larger with the least square method 
than with the step-by-step method. However, this factor 
proved to be of value in that the resulting curve was more 


realistic than the curve derived by the other method. 





Chapter VII 


SUMMARX AND CONCLUSIONS 


The research conducted here has been for the purpose of 
investigating the ຕ: “າ of utility theory as an aid for 
making business decisions. This investigation seemed valuable 
because other methods fail to incorporate the element of risk 
in a reasonable way. As it evolved, the work focused on ways 
to make the utility calculation a mechanical process. 

In order to mechanize the procedure two things were re- 
quired. One, a curve-fitting method which could be expected 
to work for all reasonable utility curves. Two, a way of 
producing data for the use of the curve-fitting process. 
These two combined would produce a mathematical expression 
for any utility curve. 

After such a mathematical function is available, an 
investigation of the usefulness of the utility method would 
be possible. Because of time limitations, this work con- 
cerned itself only with the two preliminary requirements. 

A questionnaire was designed which would generate 
probabilities from a series of questions involving two dif- 
ferent games of chance. These two games were represented 
separately as Part I entries and Part II entries. The final 
draft of this questionnaire was distributed to selected 
groups of business students at the University of Kansas. 

The purpose of the distribution was limited to the analysis 
of the effectiveness of the questionnaire as a means of 
generating useful data which could be converted to utility 


values. 
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In our case with students as subjects, the effectiveness 
of the questionnaire proved disappointing. However, those 
problem areas of misunderstanding by the students should be 
reduced considerably when the questionnaire is administered 
to business executives. With the latter subjects, the 
monetary values of the entries should not exceed the monetary 
values of the business ventures which these executives con- 
Sider on a routine basis. Likewise the interpretation of 
probability as applied to alternative business ventures should 
be understood with more ease and clarity. It was apparent, 
however, that the questionnaire is more effective when 
administered in the presence of a qualified operator. It 
can be concluded that when the problem areas of misunder- 
Standing are reduced to a minimum, the questionnaire will be 
an effective means of generating data. 

ln the development of the questionnaire it was recognized 
that individuals respond differently to different games of 
chance or types of bets. For this reason, two separate types 
of propositions were included. An effective means of in- 
corporating the data from these different propositions created 
one of the major problems in this research work. 

As noted earlier, the use of only Part I data proved to 
be the easiest. Use of Part I negative data and Part II 
positive data likewise created very little difficulty. We 
found that.participants generally placed a higher utility 
for corresponding dollar entries with the Part II proposition 
than with Part I. Because of this factor, the combined use 


of Part I and Part II data proved more difficult. 
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Two methods are presented for the incorporation of the 
Part II points with the Part I points. The first method, 
identified as the step-by-step procedure, was responsive for 
all questionnaires which were included in this portion of 
the analysis. Unfortunately, this method did not give a 
smooth curve. A modification of the method would be to re- 
calculate all previous points whenever a Part II point was 
incorporated. To recalculate each previous point by the 
step-by-step procedure, however, would be very time consuming. 
Herein lies the difficulty of this method. 

Because of the difficulties inherent with the step-by- 
step methods, we developed a procedure which would make step- 
wise calculations unnecessary. We chose to do this by 
adopting an approach whose basis was that of the general 
least square concept. A series of equations was written 
which included each of the utility values whenever they 
appeared at each data point. By solving these equations 
Simultaneously, each value of the variable was considered in 
each of the calculations. As a result there was simultaneous 
calculation of each point in the curve-fitting procedure. 
Although this method is considerably more accurate and 
Superior to the step-by-step procedure it was not responsive 
in all cases. It is believed that the inability of the 
general least square approach to respond satisfactorily to 
erratic data points caused the difficulties of this method. 
In those cases where the data points represented a relatively 
Smooth curve, this method was responsive and proved most 


Satisfactory. 
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It is concluded that our least square method is the 
more accurate and should be used in all cases where the 
program is responsive. The step-by-step procedure might be 
used in the other cases. 

The actual curve-fitting portion of the research included 
two different procedures. The general least square method 
was used to attempt to fit the curve with a logarithm func- 
tion. This procedure was found to have limited stability. 
It 1s believed that the lack of stability may be attributed 
to the same reasons that we noted earlier concerning the 
least square approach. 

An iterative procedure was used to find two values for 
the unknown variable so that the first derivative of one was 
greater than zero while the other was less than zero. This 
method, identified as the twenty question method, proved 
more useful. One possible source of difficulty is shown 
when the largest negative dollar value is quite small. In 
these cases, the iterative process of finding a value which 
will make the first derivative less than zero is quite 
lengthy. Minor modifications to the program should reduce 
this difficulty and can be made for those specific cases. 

The results of the research of the data which was 
generated by the questionnaire have not been as conclusive 
as had been hoped when the project was commenced. Although 
the questionnaire approach has avoided any direct mention 
of human behavior, this element has been thoroughly apparent 


in the probabilities. The inclusion of this factor, human 
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behavior, in the probabilities is shown both in the erratic 
progression of the data points and in the Part I and Part II 
first quadrant points. Each of these results has increased 
the difficulties of curve-fitting and obtaining consistent 
results. Because of these difficulties, conclusive statements 
are not justified, without further research, as to the type 

of curve or type of program should be the most effective 

with the questionnaire. 

Although this work has concerned itself only with the 
development of the mathematical expression for a utility 
curve, the limited analysis of the data justifies the fol- 
lowing comments. The several computer programs can be used 
to calculate the utility function; then the best curve can 
be used to plot the function. As an alternative, the actual 
utile-dollar points can be readily plotted by a Library 
computer program or easily done manually. This plotted 
utility curve has several fold advantages for the decision | 
maker. For example, a method of introducing risk as a 
substitute for pure hunch and intuitive reasoning to the 
decision is available with the use of the utility curve. 

The questionnaire-computer package offers a quick means for 
revising the criteria for decision making to keep step with 
tne rapidly changing economic status of the company or of 
the community. It offers a ready OS for subordinates 
to make decisions for management which would be consistent 
with general management policies. Likewise, management, 
itself, has a ready reference which could be used to add 


consistency to the decision. 
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This work may form the basis for further research which 
could develop methods of using the utility theory to make 
simple business decisions by computers. For further research 
in the investigation of the usefulness of the utility method, 
these recommendations are made. More extensive tests with 
the suggested least square method of data smoothing of Part 
I points and Part II points are recommended. Likewise 
extended research into the adaptability of investment type 
questions into the questionnaire is recommended. 

The questionnaire-computer package has not been perfected 
to the extent that it is the ready answer for business 
executives faced with decision making. However, the uses 
of this package as a workable aid to decision making should 
prove of assistance in delegating authority to subordinates 


and in obtaining consistency in decisions rendered. 
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Appendix I 


The purpose of the thesis is to investigate the posibil- 
ities of the use of the utility theory in decision making. 
For this reason we would like to stress that the date must 
represent the bets you would be willing to make based upon 
your personal financial position TODAY. Although we will not 
attempt to solicit an actual bet from you based upon your 
answer, this might well be the best criterion upon which to 
base your answers. 


Since the first draft, several personal factors of the 
volunteer have been deemed important to the final analysis of 
the data. It is, therefore, requested that you answer the 
additional questions listed below. This information will be 
considered confidential; any data that might be published 
would be by groups, such as married juniors, single seniors, 
etc. 


Year group in school: Marital status; Number of 
Junior Single Dependents 
Senior Married 
MBA lst year 
MBA 2nd year Male 
Other | Female 


Please place an X beside the type of housing you occupy at 
the University: Dormitory 


Fraternity or Sorority Live in property which you own __ 


Live in rental property | Live at home with parents 


Financial Status: (probably the most important question on 
the list) 

Mpproximate yearly income (include spouse's income) 

Are you paying all of your own school expenses? 

Are you on a scholarship? If so, what percent of 
your .school expenses does this represent? 

Are you receiving other financial aid while in school? 
If so, what percent of your school expenses does this TT 
represent? 

Áre you working pa part time while Е school? 


The author welcomes any constructive comments or 
criticism on this questionnaire. Such comments can be made 
on the back of this sheet. 








7? 
UTILITY CURVE QUESTIONNAIRE 


This questionnaire is part of an investigation of the 
utility theory as applied to decision making. It consists of 
a series of questions about investments. To simplify the 
Situation, however, a game of chance will be substituted for 
investments. In the questions you will be asked to estimate 
the percent probability, that your opponent will not draw a 
winning card, which you would require in order to ACTUALLY 
make the series of bets. It is hoped that these probabilities 
can be used to evaluate an investment of equal monetary value. 


To insure the greatest validity of the questionnaire, it 
is requested that you base your answers on these two criteria. 
First, your answer should reflect the EXACT probability which 
you would require to make the bet. Secondly, in the spirit 
of the competitiveness of the business world, your answer 
Should reflect the MINIMUM probability you would require. 

This latter criterion is associated with the practice of the 
lowest bidder being accepted in a call for sealed bids. 


For Part I of the questionnaire, the game will be played 
with a deck of 100 cards consisting of black cards and red 
cards. The results of the game depend upon one and only one 
card being drawn by your opponent. Should your opponent draw 
a BLACK card, you WIN; should he draw a RED card, you LOSE. 
You can deSignate the number of black cards in the total 
deck of 100 to establish the conditions under which you would 
be willing to bet that your opponent will not draw a red card 
and win the bet. It might prove helpful in thoSe cases where 
you are in doubt as to a satisfactory answer to follow the 
below listed series of questions, or Similar questions, until 
an acceptable ratio is found. 


Would you bet that your opponent would not: 


(1) draw a black card out of an ordinary deck of playing 
cards? If so, the probability in our game would be 
50, 


(2) araw a spade? If so, the probability would be 75. 
(3) draw an ace? If so, the probability would be 92. 


(4) draw the Ace of Spades? If so, the probability 
would be 98, 


In some cases you may not want to make the bet even with 
99% probability that your opponent would not draw a wimning 
card, In those cases, please indicate NO BET as your answer. 





78 
PART I OF QUESTIONNAIRE 





Поло ОГ ກວາ. Percent 
Win Loss Probabild$3 


2 |. ^ 10 
290, 100. 
000. 1,000. 
O. 103 
000. DOE 
20,000. 20.0002 
20. ГО. 
200. 0. 
4,000. 1,000. 
0. 20. 

ОО ооо. 20 OOO 
40. 40 . 
4,000. 2,0008 
ШІ 000. 70.080. 
100. 200 
BOO. 100. 
2,000. OUS 
200, 2008 
1,000. 1.00 
500. 00. 
0,000. 000. 





79 
PART IL OF QUESTIONNAIRE 


On this second part of the questionnaire ССИ СЕЕ 
ists of a series of pets on ode и 
will havo a 50/50 probability of winnie о ват үл 
indicato in the second column the amount you would expect to 
win in order to play the same on the basis of the amount of 
loss indicated in column (1). You are betting on an even 
orobability that your opponent will not win the toss on one 
and only one toss of the coin. Remember, you nave an even 
enance of winning or losing. Your answer may be recorded as 
100 to indicate a $100 bet. 


Sunt of Amount of Percent 
LOSS Win PXLODO DIA 

EM i0. = 50/100 

4,000. 50/100 

20. 50/100 

К 000. 50/100 

20,000. 50/100 

ШОО 50/100 

1 000. 50/109 

50. 50/100 

200. 50/100 

Ш Ооо. 50/100 


400. DOV ACO 
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FORTRAN PROGRAM FOR LEAST SQUARE METHOD 
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Appendix II 
COMPUTER PROGRAM TO CALCULATE UTILITY VALUES 
BY THE LEAST SQUARE METHOD 


Purpose 


The purpose of this program is to incorporate the Part 
II data points with the Part I data points. The data 
smoothing technique used is that of solving a series of 
equations simultaneously. The utile-dollar data points 
derived by this program can be used as input data for one 


of the regression analysis programs. 


Language 
Fortran IV (IBM 7040 Computer). 


Symbolic Dictionary 


VARIABLE S/A* I/O** DESCRIPTION 
N S I Total number of Part I data points. 
ка 5 I Total number of Part II data points. 
ER A I&O Probabilities used to determine а 


positive utility value from a = 
tive utility value. 


PRN A I&O Probabilities used to determine a 
negative utility value from a known 
positive utility value. 


21 A I Dollar values for the Part II data 
points. 
X A I Dollar values from questionnaire 


entries which become the dollar 
values for the positive Part I 
points. 





#5 - Single variable; A - Array of variables. 


**T - Input; O = Output. 
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U A 0 Utility values for the positive 
Part I data points. 
UN A O Utility values for the negative 
Part I£ data points. 
UNXL A = Utility value for the Part II data 
points. | 


Program Routine 


This program uses the general least square concept to 
incorporate data points from two sources. It is divided into 
two sections, namely, the matrix generater and the matrix 
solver. 

To develop the matrix generater, the following steps 
were used. А Series of equations were written which invluded 
each of the utility values (Us, UN 4, UNX1;) whenever they 
appear at each data point. The partial derivatives of each 
equation with respect to the three variables were established. 
A square array with dimensions equal to the number of equa- 
tions was used to represent the matrix generater. The rows 
within the array represented the derivatives of the indif- 
ference equations starting with e4 through ey. The columns 
represented the various variables and were designated from 
left to risht as U, through Ups UN} through УМ), UNX1, through 
UNX1L,, and the last column contained the known value U. 

(For this explanation only, n equals the nunber of Part I 
data points while m equals the number of Part IL data points.) 
Whenever & partial derivative was not equal to zero, its 
coefficients were computed into the proper array position. 
Once all of these values were entered, a library program was 


used to solve the matrix. 
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By solving these equations simultaneously, TE value of 
the variables is considered in the calculation thus there is 
Simultaneous calculation of each point in the curve-fitting 
process. 

As this program is written, all dollar values must be 
transferred to one array and ordered from largest negative 
to largest positive before this data can be read into either 
of the regression analysis programs. The corresponding 


utility values are ordered accordingly. 





MAINLINE FLOW DIAORAM FOR 
LEAST SQUARE PROGRAM 


50 
Ie 1, 
J 


K)PR(K) + PRN(K) ` 


p^ PR(K)PR(K) 
+ 
) = (1 - PRN(K))PRN(K) 
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200 
K*1,N 









KK = K +H < 1 

KI = K #1 

KKK = Of] * K <ü 1 

A(KK,K)=(PR(KI))e(1, = PR(KI)) 
N = N < 


5 


Л (КК, КК )<А (КК,КК) + (1, - РА(КГ)) «2 





T 


F 


A(KK,KK)=A(KK,KK) + 1, 
A(KK,KK)=A(KK,KK) + (1, - PRN(K))ue2 







A(N,LI,)==(PhN(1))oe 
(1. = PHX(1)) 








KK . К + ໂ = 1 
KKK = 2* N+ K = 1 
A(KK,KKK) = 1 





9 


LJ = LLL(J 
Ll = LJ - 1 
K = LJ + 1 
















DELX e X(LJ + 1) - X(LJ) 
DX = (X(K) - x1(J))/DELX 
DXX = (Xi(J) - X(LJ))/DELX 
ММЈ = 2*N + J - 1 



















rl D) 4 DX**2 
L1,1J) = A(L1,LJ) + DX*DXX 
at - rd - DX 
NNJ,L1) = A(NNJ,L1) - DX 
LJ,L1) = A(LJ,L1) + DX*DXX 





> > > > > 









A(1,LL) = A(1,LL) 

- DXX*DX 
A(NNJ,LL) - 
A(NNJ,LL) + DX 


nen = A(NNJ,LJ) - DXX 
A(NNJ,NNJ) = A(NNJ,NNJ) + 1 


E] 





5 


KKK = N + K - 1 
UN(K) = COE(KKK) 


oy 


600 
K=1, NX1 t 


WRITE 
U(K), UN(K), UNX1(K) 


GO TO 1 
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FLOW DIAGRAM FOR SUBROUTINE EMSLVR 


AX = A(KKK,I) 
А(ккк,1) = 


A(IBIO,I) 
AC IBIG, I) AX 










AY = A(J,KKK) 
ا‎ 
A(J,JBIG) = АХ 


A(J, NP1)e -A(J,NP1) 








IDUM e IC(KXX) 
рК « IC(JBIO) 
IC(JBIO) - IDUM 





0 = 


CONTINUE 


NN=NN + 1 
NNN=NN + 1 


ммемм + 1 





CHECK FOR SINGULAR MATRIX 


00 
I=NN, N 


PRINT SINGULAR MATRIX 


1900 


' = 1,N 
<> ACOE(I) =O. 
z 
Bo н 
J=NNN, NP1 


1900 
А(1,2) = A(1,J3)/A(J,NN) 


CONTINUE 






AB = ABS(A(1,3)) 
IBIQ = I 
JBIG = J 


CONTINUE 


0 e 
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I=NNN, 


SUM = SUM - A(J,I)COE(I) 








БФ———— д 


К = x 
ACOE(K1) # COE(I) 


CONTINUE 


CONTINUE 


WRITE 
VARIABLE IDENT 


WRITE ACOE(J),J,IC(J) 


U 
g 
sus 
4 
' 
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10 


20 
25 


50 


= 


40 


41 


30 
200 


DO 200 К + 1 N 


88 


THE LEAST SQUARE PROGRAM DESIGNED TO INCORPORATE 
PART II POINTS BY USE OF SIMULTANEOUS EQUATIONS, 


ME - = — 


DOUBLE PRECISION A»sPR»PRN»DELX»X»DX»DXX»UsCOE»UN»UNX1»X1»P 
DIMENSION UNX1(50) » СОЕ! 5О ) +) ЕГ (50) »Х1 (50) »УХ (50) 
DIMENSION PR(50)» PRN(50)» P(50)» A(50551) » 50)» ИМ5О) 
READ (594) N 

READ (594) МХ] 

READ (593) (Х(Т)» I -1»М) 

READ (593) (Х1(Г)» I =1»NX]) 

READ (593) (PR(I)>» I = l» N) 

READ (593) (PRN(I)» I = l»sN) 

READ (5.6) №1» №» М3» Ма» М5» М6 

FORMAT (6A6) 


AAA 


U(1) = le 

DO 10 J = 1» МХ1 

LLL(J) = LOOK (X»X19»J) 
IF(LLLUJ)eGE.N) GO TO 20 
CONT INUE 

GO TO 25 


RAL = J = 1 


KI = K 
А(КІЭКІ 
KK = K + N 
А(КІЭКК) = 
KKK = KK + 
А(КІЭККК) 


)*PR(K) + PRN(K)*PRN(K) 
- РК(К)) “РК (К) 
(le — PRN(K))*PRN(K) 


KK = K +N - 1 

KI = K + 1 

A(KKsK) - (РК(КІ))%(1» - РК(КІ)!) 

ККК = 2*N + ເ - 1 

KN = N - 1 

IF (KeGTeKN) GO TO 40 

А(ККЭКК) + А(ККУКК) + (le - РК(КГ)) #2 
IF(K«GT.NX1) GO TO 41 

AC(KKsKK) = A(KKsKK) + le 

AUKKSKK) = A(KKsKK) + (lo - РЕМ(К)) #2 
IF (KeLTe2) GO TO 30 

КККК = К = 1 

A(KKsKKKK) = PRN(K)* (le - PRN(K)) 

GO TO 200 

A(NsLL) =-(PRN(1))*(lo - PRN(1)) 
CONTINUE 
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300 


&50 


149 


| 
297 
29% 


700 


34 


800 


35 


1105 
33 
449 
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DO 300 K = 1, NX1 


KK = 2 * N + K - 1 
A(KKsKK) = le 

KKK = N+ K - 1 
А(ККЭККК) = le 

DO 450 K = 19NX1 


KK = K + N - 1 
KKK = 2*N + K - 1 


A(KKsKKK) = 1, 
DO 299 4 51» МХ1 
LJ = LLL(J) 

pr te I 

S E Ad El 


IF (LJeGTe6eN) GO TO 333 
DELX = X(LJ + 1) - X(LJ) 


DX = (X(K) - X1(J))/DELX 


DXX = (X1(J) - X(LJ))/DELX 
ММУ = 2#N + J - 1 

IF о ЕО 1 ) GO TO 135 
A(L1sL1) = AlLlsL1) * DX*x2 
A(L1sLJ) = A(L1sLJ) + DXF*DXX 
A(L1sNNJ) = A(L1sNNJ) - ОХ 


А(ММ./2(1) А(ММ./9(1) - DX 
A(LJsL1) = A(LJsL1) + DX*DXX 
A(LJsLJ) = A(LJsLJ) + ОХХЖЖО 
ACLJsNNJ) - A(LJSsNNJ) - DXX 
AK(NNJSLJ) 7» A(NNJSLJ) - DXX 
А ( ММ.) ММ.7)) < А( ММ.) „ММУ) + le 
GO TO 299 

A(1sLL) = A(1sLL) -DXX*DX 
A(NNJsLL) = A(NNJsLL) + DX 
GO ТО 149 

CONTINUE 

DO 700 1 = loLL 

WRITE (65900) (А(І».))» J = 1912) 
CONTINUE 


WRITE (652) 
FORMAT (1H1) 
IF (LL - 13) 33933934 


DO 800 I = lett 

WRITE (63900 ) (AC TsJ)9 J = 13924) 
CONT INUE 

WRITE (692) 

IF (LL - 25) 33. 33» 35 

DO 1105 I = 1sLL 

WRITE (6900) (А(1»))» J = 25936) 


CONTINUE 
NP1 = 2*N + NX1 - 1 
CALL EMSLVR (AsCOEs1sNP151) 
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400 


500 


600 


1100 
900 


10 
20 


90 


00 400 К = 2. М 

KK =К - 1 

U(K) = COE (KK) 

DO 500 K = l9 N 

KKK = N + K - 1 

UN(K) = COEIKKK) 

DO 600 K = 1sNX1 

L = 2 *N + K - 1 

UNX1(K) = COE (L) 

WRITE (7s8) №]. М2» М3» М6 


FORMAT (36Xs3A6»20X»A6) 

WRITE (691100) (0(%)» K = 19N) 
WRITE (691100) (UN(K)s K = 19N) 
WRITE (6,1100) (UNX1(K) + K = 1sNX1) 
WRITE (791100) (U(K)»s K = 19N) 
WRITE (751100) (UN(K)s K = 13N) 
WRITE (791100) (UNX1(K) +s K = 1sNX1) 


FORMAT (1Х» 6(Ғ8.2»2Х)) 
FORMAT (1Х» 12Ғ10.7) 


GO TO 1 

FORMAT (F10.3) 
FORMAT (12) 
END 


DOUBLE PRECISION FUNCTION LOOK (XsX1»J) 
DOUBLE PRECISION X»X1 

DIMENSION X(50) X1(50) 

DO 10 L = 1,100 

IF (X(L).GT.X1(3)) GO ТО 20 

CONTINUE 

LOOK = - 1 

RETURN 

END 
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SUBROUTINE EMSLVR (AsACOE sMORDER 9NoMPTO) 


11 ARDER LLIE: MATRIX BEEORE .ELACH.UELIMINATION:LE MORDER=+1 L — 
N= ORDER OF MATRIX 
WILL SOLVE AN (N)X(N+1) MATRIX 
REQUIRES MATRICES OF THE FORM (A)X(COE)=(B) 
ACOE=VARIABLES TO BE SOLVED FOR 
A(IsJ)= MATRIX ENTRIES 
COLUMN (I»sN*1) OF A MATRIX CORRESPONDS TO- COLUMN- MATRIX B . 
DIMENSIONED VARIABLES MUST BE AT LEAST OF ORDER N OR N+1 
DIMENSION A(NsN+1)»s IC(N)»s COEIN+1)s АСОЕ М 
ANSWERS TO SINGULAR MATRICES ARE ZERO(O) 
DOUBLE PRECISION ACOE sAsCOE  sAXsAYsSUM 
DIMENSI BE A C OE( (50), A505 IC 50s COE (5L) 
NM=N = ວ. 
NN =0 
KK =0 
MM =0 
NP1=N+1 
NM1=N-1 
DO 3 J=loN гы E ະຫນ 
3 AL JsNP1)z-A(JS NP1) 
“ INITIALIZE SUBSCRIPT COLUMN 
799 DO 800 ,/:1»М 
ຍວ dl ER 
KKK=0 
IF (MORDER) 9995755999 ل‎ ы وک‎ 
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999 KKK=KKK+]1 
АВ <ОАВО (А (ККК ,ККК ) ) 
IBIG=KKK x ЩЕ ——— — 222 —— 
JBIG=KKK 
DO 901 I=KKKoN 
DO 901 J=KKKoN 
IF (AB-DABS(A(1,J3)))90099019901 

900 AB =DABSI(A(I>J)) 
IB I6=1 
JBIG=J 

901 CONTINUE 

910 DO 920 I=1» NP] 
АХ=А(ККК»1) 
А(ККК»І)-А(ІБВІб»іІ) 

920 A(IBIGsI)=AX 

DO 930 J=1»N 
AY=AUJ»KKK) 
AULJ»KKK)=SA(J»JBIG) 

930 AlJ»JBIG)=AY 
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CY CY Y 


OOO 
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CY CO OY 
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940 


4 
Ше 


ner 
79 


80 


85 


90 


91 
92 


100 


109 
110 


IDUM=IC(KKK) 
IC(KKK)=IC(JBIG) 

IC (JBIG)=IDUM 

IF (NM1-KKK) 715715999 
CONT INUE 

ММ =ММ+1 

NNN=NN+1 


 MM=MM+1 


чо» ало» чо» чи» «и» ແນ ແບ ແພບ ap чо» op s G s G “s G ш» Фф ow s G G «ні» «ы» тəн» «ні» е» еки» ы» G G (uum amb 


Ld JE UL سے سے جت سے ج‎ ТЕТ Я 


IF (A(NNsNN)) 7791700977 


DO 80 I=NNoN 

IF (ACI sNN)) 795980979 
DO 80 J=NNN>NPI 
ACISJ)SACISJ)/ACI 9 NN) 
CONT INUE 

KK=KK+1 

IF (KK—-NM1)85 9859100 
DO 90 I=NNNoN 
IF(A(TINN)) 89990989 
DO 90 J=NNN>NPI 


A(T9J)=A(19J)-A(NN9J) 


CONTINUE 

TF (MORDER) 91975991 
IF (KK=NM14+1)92992975 
KKK=MM 

GO TO 999 


чо» чт» G. G. G. s чи» ow GERD ແຫນ ແແ а» G q ча» чо» G Фи» чо» чо» s чо» что G G G s ap Ab чо» = a= == 


CHER cum ime “< Q сш» J W s GNEND m q“ СЯР «ММ Qm “WE сольно скс, mam Qa Q q ss Q m O um s ae qa G m s dini 


СОЕ (МРЪ) 5140 

ОО 110 К-1»ММ 
50М-0.0 

J=NP1-K 

L=J+1 

DO 109 I=L»NP] 
SUMzZSUM-A(JsS I)*COE(I) 
COE (J)=SUM 


-—-— .ແແບ ແຫນ ແຫນ ແບບ G s Q ແຫບ ແນບ ແຫບ ແຫບ 3 Q ຈາ ແຫນ ແຫນ өнер ແຫນ ແນບ ແນນ ແນ ແນ s“ «жы» ແນ A омы» екн» «мн» ແເ Q 


DO 1005 151 „ММ 
K1=1C(1) 
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| зу} Тм Сао 






IE AAA 
E Tecum 
IRE 9903 тин" | 
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1005 


1500 


1599 
1600 


1700 
10 


1900 


ACOE (K1)=COE (I) 
CONTINUE 
ІР (МРТО ) 150091600 • 1500 
CONTINUE 
WRITE (692) 
DO 1599 Ј=1,ММ 
WRITE (691) ACOE(J)»* Js ТС(.)) 
CONTINUE 
RETURN 
ЕОКМАТ(1Х»ІН Е15.6>218) 
FORMAT(1Xs1LHO41H VARIABLE IDENT 
PRINT 10 
FORMAT (1XslLHOs1l6H SINGULAR MATRIX) 
ОО 1900 Із1»М 
АСОЕ (1) 506 
RETURN 
END 


oo 
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APPENDIX IIE 


FORTRAN PROGRAM FOR TWENTY QUESTION METHOD 


94 
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Appendix III 


COMPUTER PROGRAH TO CALCULATE 
COEFFICIENTS FOR A UTILITY CURVE 


Purpose 


The purpose of this program іс to calculate the coef- 
ficients for a utility curve from the utile-dollar data 


points computed by one of the data smoothing programs. 


Language 
Fortran IV (IBM 7040 Computer). 


Symbolic Dictionary 
VARIABLE — S/A* . I/O** | DESCRIPTION 


N S I Total number of data points. 

X A I&O Dollar values for the data points. 
U A I&O Utility values for the data points. 
B 5 Ó The unknown variable within the 


logarithm function which causes the 
curve to go through the (0,0) point. 
Initially it is a guess to start 
processing. Thereafter it is 
computed internally. 


BG 5 - Tnat value of the variable B whose 
first derivative is greater than 
Zero. 

BL 5 - That value of the variable B whose 


first derivative is less than zero. 


BM 5. - That value of the variable B whose 
first derivative is equal to zero or 
within limits set by the progran. 





#5 = Single variable; A =- Array of variables. 
##I — Input; O - Output. 
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UCAL A О Calculated utility value for tne 
data points. Computed internally 
with the use of the computed 


coefficients. 

DEVIAT A е; Deviation of data point utility and 
the calculated utility values. 

STDEV S О Standard deviation of the data point 
utility and the calculated utility 
values. 

A S Ó The unknown variable which controls 


the magnitude of the curve. This 
is computed internally. 


MN S I Number of iterations in the process 
of calculating BG's and BL's to 
determine BH. 

MTEST S I Total number of X values dropped. 


NBTEST 5 I Number of iterations to find an 
FOFB less than zero. 


Program Routine 


In this method of iteration, two values of B are deter- 
mined so that the root of the equation lies between then. 
Since the root lies between the two numbers, the graph of 
FOFB = f(B) must cross the x-axis between B - B4 and B = Bo, 
and the FOFB, and FOFB, must have opposite signs. We designate 
the value of B which gives FOFB less than zero as BL, while 
the value of B which gives FOFB greater than zero as BG. 
Having established these two limiting points, BL and BG, we 
used 20 iterations to calculate the correct value of B 
(designated BM) where: 

BM = (BG + BL)/2. 


The method of determining the FOFB follows the general 
mathematical model described in Chapter III. In brief, the 





27 


general mathematical equation was set equal to е. Botn 

Sides of the equation were squared and the partial derivatives 
of E with respect to a and b were set equal to zero in order 
to minimize the total error. Several collective terms which 
were of a repetitive ae Ha designated by alphabetic 
titles and included in the program as function programs, This 
allowed FOFB to be identifled in the malnline as & mathemati- 
cal expression consisting of the several functions. 

When BM is determined, this value is used to calculate 
the other unknown variable A. The program calculates the 
utility values for the dollar points using the calculated 
coefficlents. These values are used to determine the standard 


deviation for the calculated utility curve. 











FLOW DIAORAM FOR 


TWENTY QUESTION PROGRAM 









ORDERED DOLLAR AND 
UTILITY VALUES CAN 
BE READ IN FROM 
EITHER OF THE DATA 
SMOOTHING PROGRAMS 








J» j +} 
CALL CLOCK (J) 





TRANSFER X(I) TO X1(1) 


455 
I = 1,N 






SAM(J) = ALOG(XJB) 


FOFB1 = T(X,N,B,U)» 
W(N, SAM) = Q(N, SAM, 
U, )*R(X,N,B, SAM : 





> + - 


BTRY e - x(1)*2 


CONTINUE 





SAM(J)=ALOG(XJB) 


FOFBT e T(X,N,BTRY,U 
W(N,SAM) - Q(N,SAM,U 
R(X, N, BTRY, SAM 


WRITE 
BT, FOFBT 








) 





XJB-((X(J) + BTRY)/BTRY) 










FOFB(1) = FOFB1 


KK = KK + 1 
В = В ~ 150 














FOFB(1) = FOFB1 
PTEMP = B 

JJ = JJ + 1 
В = В = 150 


OO TO 19 


FOFB(II).aT.0.0~e<Porn(1).07.0.0 


Р 
BG = PTEMP 
T BL = B я 


OFB(I).aT.0.0 





00 то 905 
В É 
DL e NTEMP O 
T 
FOFB(II).QT. FOFB(I]» 
GO TO 905 





BTRY » BIRY - X(1) 


NDTRY = NBTRY + 1 B e NTEMP - 150 


(9 LP KK = KK + 1 
B - 15 


Ho 


B = 1 
KK = KK + 


t 


; 


1 


X 






WRITE 
NO FOFD FOUND GREATER THAN ZERO 





FOFB(II). LE. FOFB(I) 







B = PTEMP + 100 
JJ = JJ + 1 


PTEMPT e B 
B = В + 100 
JJ = JJ + 1 


<> т <> 
ре 


| < 






RE RED 
IT = I+1 BO e PTEMP 
(00) BL = NTEMP 
FOFB(I) = T(X,N,B,U)* (к 1.) NCHECK = O 





PUN) - Q(N, SAM,U)* 
R(X,N,B, SAM) 











=“ K + 1 


x € TRANSFER XA (I) TO X(I) 
- па 


x(n) 








° @ 
E. 


1777 г 


I =i, 
ສສກ > 














MCHECK e MCHECK + 1 






WRITE 
CALCULATION DID 
NOT CONVEROE 





WRITE 
1, BM, FOFBM 


NCHECK = NCHECK + 1 


i - 
F 





FOFBM. LE. 0. O 


© 
WRITE 
XDATA, UDATA, UCALC, DEVIATION 


| 


UCALC(J) = A(ALOO(X(J) + B)/B) 


DEVIAT(J) = U(J) - UCAL(J) 





j WRITE 
X(J), U(J), UCAL(J), DEVIAT(J) 


600 


. 
| (1-2) 
1 


SUM = SUM + (DEVIAT(I))(DEVIAT(I)) 





100 





ам 


433 
18 


2300 


E 


13 


2400 


12 


3300 
3301 


3302 
2596 


3400 


101 


THIS SECTION OF THE PROGRAM IS THE TWENTY QUESTION METHOD 
OF REGRESSION ANALYSIS. IT FOLLOWS THE ORDERING SECTION, 


МСНЕСК + 0 
NBTRY = O 


TRANFER X(I) TO X1(1) 

DO 433 I = 1sN 

X1(1) = X(1) 

В = - Х(1) % 100» 

DO 2300 J = 1»М 

XJB = ((X(J) + B)/B) 

SAM(J) = ALOG (XJB) 

FOFB1 = T(X+sNsB+sU)%W(NsSAM) — Q(N+SAM+U)*R(X+N+B +SAM) 
WRITE (6511) B» FOFB1 

FORMAT (1X»4H B =>F10.3>10X>7HFOFB1 =>E15+6) 

BTRY = 15000, 

CONTINUE 

DO 2400 J = 19N 

XJB =-((X(J) + BTRY)/BTRY) 

SAM(J) = ALOG (XJB) 

FOFBT = T(X*NsBTRYsU)%W(N+sSAM) - QIN»SAM»U)J#R(XsN»BTRY»SAM) 
WRITE (6512) BTRYsSFOFBT 

FORMAT (4H BT=>»F10+.3»10X>»7HFOFBT =>E15+6) 

IF (FOFB1+GT+0+0) GO TO 3400 

IF (FOFBTeLEe0.0) GO TO 3300 SS — 


BG = BTRY 
BL = B 
GO TO 905 


IF (FOFB1eGTeFOFBT) GO TO 3600 

BTRY = BTRY - Х(1) 

NBTRY = NBTRY + 1 я 

IF (ABS(FOFBT)+.LE..00001) GO TO 3302 
IF (ABS(FOFB1)+.LE+.+.00001) GO TO 3303 
IF (NBTEST - МВТКҮ) 15» 15» 13 

WRITE (699) 

FORMAT (32H NO FOFB FOUND GREATER THAN ZERO) 
OU FPO T 

BM = BTRY 

GO TO 1800 

BM = B 

GO TO 1800 

IF (ЕОЕВТ» Г Е»0,0) GO TO 3500 

IF (FOFB1eGTeFOFBT) GO TO 3301 
FOFB(1) = FOFB1 

PTEMP = B 

B = B — 150,» 

JI" > JIJ +1 
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GO TO 19 
3500 BG = B 
Bi TR 
GO TO 905 
3600 FOFB(1) = FOFB1 
NTEMP = B 
B = B — 150, 
KK = KK + 1 


GO TO 19 
"19 Nr 2 
20 CONTINUE 
DO ft ie = ІТ» М 
DO 300 J s l»N 
XJB = ((X(J) + B)/B) 
300 SAM(J) = ALOG(XJB) 
FOFB(I) * TUXGNSB»U)*W(NSSAM) — Q(NsSAMoU) #R(X9N9B 9SAM) 


WRITE (6916) I»B»FOFB(I) 
16 ҒОЕМАТ(ЗН 1-»12»10Х»З3НВ =»>»F10.3»10X»6HFOFB -»Е15.6) 
II  I- 1 
IF (FOFB(II)+.GT.0.0) GO TO 720 
IF (FOFBII).«4GT.,0O,0) GO TO 200 
TF (FOFB(II)-eGTeFOFB(I)) GO TO 710 
B = B - 1509 
KK = KK + 1 
GO TO 776 
710 B = МТЕМР - 150 
кк - КК + 1 
GO TO 776 
200 BG = B 
BL = NTEMP 
GO TO 905 
720 IF (FOFB(I).GT.0.0) GO TO 201 
BL = 5 
BG = PTEMP 
GO TO 905 
201 1Е (РОРВ (ТТ) ЕЕ, РОРВ(Т))} СО ТО 203 
РТЕМР = В 
В = В + 100. 
JJ, z JJ +<) 
GO TO 776 
203 B = PTEMP + 100, 
JJ JJ * 1 
776 IF (KKeGEel) GO TO 820 
778 IF (B + X(1)) 810, 810, 777 
GO TO 777 
820 IF (JJeGEel) GO TO 920 
GO TO 778 
777 CONTINUE 
780 WRITE (65790) 
790 FORMAT(30H CALCULATIONS DID NOT CONVERGE) 
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CY CÓ OY 


810 


17 


920 


905 


334 


910 


1300 


14 


999 


1000 
1777 
1800 

500 
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GO TO 1 
N = N-1 
IT = I + 
рО K 
KKK = K + 1 

X(K) = X(KKK) 

U(K) = U(KKK) 

MCHECK = MCHECK + 1 

IF (MCHECKeGTeMTEST) GO TO 780 
GO TO 20 

BG = PTEMP 

BL = NTEMP 

NCHECK = O 


TRANSFER X1(I) ТО Х(Г) 


«к «ам» ED ແຫນ ແຫນ «бы» «ва qub qub CED CD Guam Фа» «па» ФИШ» Сс bé enm ແຫນ “QE m ¿Qm ແທ 


DO 334 I = 1+sMM 


1 
= 1»9М 


X(1) = X1(1) 
DO 1777 I = l1» M 
BM = (BG + BL)/2e 


DO 1300 J = 1» ММ 

XJB = ((X(J) + BM) /BM) 

SAM(J) = ALOG (XJB) 

FOFBM - Т(Х»М»ВМ»/)) «М(М»5АМ) - 0(М»5АМ»/) ЖК(Х»МэВМ»5АМ) 
МКІТЕ (6914) І» ВМ» ҒОЕВМ 

FORMAT (3H 1=s12»10Xs4HBM =sF10.3»10X,» THFOFBM =sE15.6) 
IF (ABS(FOFBM) eLE«e000001) GO TO 1800 

NCHECK = NCHECK + 1 

IF (NCHECKeGTeNTEST) GO TO 1800 

IF (FOFBMeLE.0.0) GO TO 1000 


BG = BM 

GO TO 1777 

BL = BM 

CONTINUE w 
B = BM 

A = T(XsNsB+sU) / К(Х»М»В»5АМ) 


Al = Q(NsSAM+sU) / М(М»5АМ) 
WRITE ( 69700) 
WRITE (69800) As Als B 


. WRITE (65,151) М1» М2» N3 


151 


600 


30 


FORMAT (/36X93A6/) 

WRITE ( 69850) 

DO 600 J = 1» MM 

UCAL (J) = A¥(ALOG((X(J) +B)/B)) 


DEVIAT(J) = U(J) - UCAL (J) ga 
WRITE (65900) X(J)» UtJ)» UCAL(J)» DEVIAT (J) TTT 
SUM = 0,0 


DO 30 I = 1» ММ 
SUM = SUM + (ОЕУТАТ (Г) ) 2 (ОЕМТАТ (ГТ) ) 
AM = MM - 1 
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STDEV = SQRT (SUM/AM) 
WRITE (6931) STDEV 
31 FORMAT(22H STANDARD DEVIATION =9F10-5) 
700 FORMAT (1H1931X918H UTILITY CURVE FIT//) EA | 
800 FORMAT (7Xs3H A=-»F10.3»10X»4H Al=-»F10.3»10X»s3H B=-,F10.3) 
850 FORMAT (7Х» 5HXDATA* 15X$55HUDATAs 15X$55HUCALC» 12% 
19HDEVIATION) 
GO TO 1 
WRITE (692) 
2 FORMAT (1H1) 
900 FORMAT (4X» 3(F10.3»10X)» E12.6) 
END 


FUNCTION @ ( N» SAM» U) 
DIMENSION SAM (20)» U (20) 


Q = 0.0 
DO 10 J= 1» N 
10 Q =Q + U(J) * SAM (J) 
RETURN 
END 4 


FUNCTION R ( X » М» В» ЗАМ) 
DIMENSION X (20) SAM (20) 


R = 0.0 
DO 10 J= 1» N 

10 R = R + SAM(3Ə Y * ( 1lə / (X(J) + B)-1e/B) 
RETURN 
END 


FUNCTION T ( X» N» B» U) 
DIMENSION X (20)» U (20) 
"i = 0.0 
DO 10 Jj = 1» М 
WF T'[UU Ol Le Ü J 
RETURN 
END 


FUNCTION W ( М» ЗАМ) 
DIMENSION SAM (20) 


и = 0.0 
DO 10 J = J], N 

10 М = М + SAM(J) %%2 
RETURN 


END 
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APPENDIX IV 


FORTRAN PROGRAM FOR STEP-BY-STEP METHOD 
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Appendix IV 


COMPUTER PROGRAY TO CALCULATE UTILITY VALUES 


BY THE STEP-BY-STEP PROCEDURE 


Purpose 
The purpose of this progran is to incorporate the Part 
II data points with the Part I-data points by means of a 


data smoothing process. 


Language 


Fortran IV (IBM 7040 Computer). 


Symbolic Dictionary 





VARIABLE “າໍ that DESCRIPTION 
H 5 D Total number of Part I data points. 
ка 5 ШІ Total number of Part il date сорта” 
p A 10 Dollar values fron questionnaire 


entries which become the dollar 
values for the positive Part I 
ЮО ЕЗ 


PR A 120 Probabilities used to determine a 
positive utility value from a 


= 


negative utility value. 

PRN A 140 Probabilities used to determine a 
negative utility value from a 
positive utility value. 


pp d L Dollar values for the Part II data 
pointsi 


“S - Single variable; A - Array of variables. 
BEI - Inputs O - Output. 
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DX А - Dollar values for negative Part 
I points. 
U A 0 Positive utility values for 


Part I data points. 


UN A 0 Utility value for the negative , 
part I data points. 


Program Houtine 
This program utilizes probabilities derived from a utility 


questionnaire to compute utility values corresponding to the 
dollar values of the data points. The procedure used to 
calculate utility values is called the step-by-step procedure 
from the method of incorporating the Part II data. The | 
utility values are determined for Part I points using Equa- 
tion (8). After each value for Part I is calculated, the 
program tests to see if the positive dollar value is greater 
than the first Part II point. When this is found to exist, 

the normal procedure is terminated temporarily to incorporate 
the Part II point. The function is considered linear between 
the two bracketing Part I points and the Part II point. А 
mean value for the upper Part I point is calculated. This 

mean value is used to continue the process. After all calcula- 
tions, all dollar values are transferred to one array and 
ordered ascending from the largest negative value.  Correspond- 
ing utility values are ordered accordingly. 

All dollar values must be transferred to one array and 
ordered from largest negative to largest positive before this 
data can be read into either of the regression analysis pro- 
grams. The corresponding utility values are ordered accord- 


ingly. 


- е е 
K. 
-— - 
> Е > ж- = = = 
=> = = = 





FLOW DIAQRAM FOK STRP«STKP-CALCULAT IONS 









HEAD 
Ра, Рин, NI, №2 
Хі, Р, Н) 


ASSUME » 310 HAS + 1 UTILITY, 80 


и(1) е 1.0 


ШІП š 5 Ua jem 7 « PRN(1)) 





ЗТАКТ STEP-BY- STEP CALC'ILATIONS 


(ойға - PR(J))/PR(J) 
າວ e FPRN(J )) 





UCHECK(JA) = -UN(JB) 

UNEW (J) = ນ 
#P(JJ)))we(P(J)-P(JJ))*U(JJ))+U (J 

JB ® JA +} 


X1(JA! 
)/2 














UCHECK (JB) = -UN(JB) 
UNEW1(J) = ,S%(((UCHECK(JB)-U(JJ 

-P(JJ)) )e(P(J) -P(JJ )) *U(JJ) 
UNEW(J) * UNEWl1(J) 





) 
) 









ОСРЕСК (16) = -UN(JC) 
UNEW2(J) = ,S*(((UCHECK(JC)<U(JJ))/(X1(JC) 
-P(JJ)))etP(J)-P(JJ) )*U(JJ )) *uhEW1(J )/2 


UNIW(J) = UNEW2(J) 
JA = JC 


UIJ UNEN (J) 

UUIEW(JJ) # <U(J)epPR(J)/(1,=PR(J))) 
UN(JJ) ທ UNNEW (JJ) 

ункка (г) * -U J ДУА ЕРІДІ? 
(7) @ UNIER (J 

JA = JA * 1 


CONTINUE 
(>) 











WRI TE | 
Р(Т), РН(1), РХ(І), РЕМ(І) 


CONTINUE 


© 
мәхет-і 


AMS(UI1(J)).GT.3000, 


CHHTILVE 


00 ТО АМ ORDERING ROUTINE 
WHICH WILL ORDER THE. DOLLAR 
VALUES AND CORRESPONDING 
UTILITY VALUES FROM LARNEST 
NEGATIVE, INSERT ZERO VALUE, 
TO LARGEST POSITIVE VALUE, 
NOW DATA IS 1N PROPER ORDER 
TO FNTER EITHER REIKESSION 
ANALYSIS PROGRAM DESIRED, 
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CY C 


аб боба 


CY CY CY C 


CY C OY 


LI 


150 


READ (59,3) (PRN(I)»> I = 19M ) 
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THIS SECTION- OF THE PROGRAM IS-THE STEP=BY=STEP METHOD 
OF DATA SMOOTHINGe 


PROGRAM DESIGNED TO DERIVE UTILES FROM QUESTIONNAIRE 

THE P(I)S ARE DERIVED FROM PART I OF QUESTIONNAIRE 

P(1) MUST BE THE LARGEST NEGATIVE DOLLAR VALUE 

THE X1(I)S ARE DOLLAR VALUES FROM PART II OF QUESTIONNAIRE 
PR VALUES ARE PROBABILITIES FOR THE 2 TO 1 RATIO QUESTIONS 
PRN VALUES ARE PROBABILITIES OF THE EQUAL RATIO QUESTIONS 
M - NUMBER OF DATA POINTS 

DIMENSION Х1(20)» Р(20)» 5АМ(20) 

DIMENSION UCHECK (20) sUNEW(20) sUNNEW(20) sUNEW1(20) sUNEW2(20) 
DIMENSION U(20)s UN(20)s PR(20)s PRN(20) 

DIMENSION PX(20)» PTEMP(20) | = к 8 
DIMENSION UCAL(20)5» DEVIAT(20)» UUTEMP(20)» X(20) © 
READ (594) M 

READ (5,4) NX1 

READ (553) (Р(1)» І =1»М) 

КЕАО (53) (Х1(!)» І -1»МХ1) 

READ (553) ІРЕ(І)» I = 19M) 

READ (59150) М1» М2» М3» М4» М5» М6 

FORMAT (6A6) 

ASSUME +510 HAS +1 UTILITY, SO 

DETERMINE UTILITY FOR -$10 BY 

U(1) - l e0 

UN(1) =-U(1) * PRN(1) / (lo - PRN(1)) 

PX (1) -Р(1) 


“< N“sum “mam qm “s вар «вар чат» < < a a. DE ммм» «мш» «мн» «мн» ашы» онын а» чт» “mam “ms. “ma. “m ағы» «мн» «ми» «ми» чин» «ын» “man. “mms < < Чо ча» чил» Ч < яны» хыт ЧЕ» шы COND пш 


SMART WO fewer FOR (1) +X TO —2X AND (2) +2X ТО - Х 


JA = 1 
DO 101 J = 2» M 
Зекен.) — -1 
Ва н ГЕ. ЈОЈ ИНС (о р Вун ЫСО 
UN(J) = -U(J) ж РТМ(Ә) / (1s — PRN(J)) 
со Рат 9 -Р = +. ແວ ຈ“ <. #? ча ks 


12 


IF (JAeGTeNX1) GO TO 107 

IF (P(J)eLTeX1(JA)) GO TO 107 

UCHECK(JA) = -UN( JA) 
UNEW(J)=65%*(((UCHECK(JA)-U(JJ))/(X1(JA)-P(JJ)))*(P(J)- 


LEUTE UG...) +. Аша.) Д2 е 
POTII) FOTI Zee ເ. — 


WRITE (6912) JAsU( J)» UNEW( J) 
FORMAT(5HOJA =912910Xs8H U(J) =9F15e3210Xs11H UNEW(J) => 


Muse.) 
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JB = JA + 1 
IF (JBeGTeNX1) GO TO 102 
IF (P(J)esLTeX1(JB)) GO TO 102 
UCHECK(JB) - -UN(JB) 
UNEW1(J) = $.5*((CUCHECK( JB) CUCJJ) )/(X1C€ JB) -PC JJ) ) ) *(P( J) - 
ТР. 99) +III) c CONCW (27/28 
UNEW(J) = UNEW1] (J) 
WRITE (69302) JBsU(J) sUNEW (J) 
302 FORMAT (5HOJB =912910X98H U(J) =»F15.3+10X>11H UNEW(J) => 
SIT De) т. D 
JA = JB 
ат JB + 1 
IF (JCeGTeNX1) GO TO 102 
IF UP UI) seb ome (Ue) GO oho” OZ 
UCHECK (JC) = -UN( JC) 
UNEW2(J) = s5*( ((UCHECK(JC)-U(JJ))/(X1(J2C)—-P(JJ)))%(—(P Y= 
1P(JJ)) * UtJJ)) + UNEW1(J)/20 
UNEW(J) = UNEW2(J) 
WRITE (653023) JC 9U( J} sUNEW (J) 
303 FORMAT (5HOJC =912910X98H U(J) =»F15.3,10X>11H UNEW(J) =, 
1F15e3) 
JA = JC 
102 Ч(Ј) = UNEW(J) 
UNNEW(JJ) =-U(J)*(PR(J)/(1es—-PR(J))) 





WRITE (6913) JJ sUN( JJ) sUNNEW( JJ) 

13 FORMAT (5HOJJ =s12910Xs8HUN( JJ) =>F15.3>10X»11HUNNEW(JJ) => 
ISE 
UN (JJ) = UNNEW( JJ) s ae xu 
UNNEW(J) = -U(J)*PRN(CJ)/(1.« - PRN(J)) 
WRITE (6914) Js ОМО)» UNNEW(J) 

14 FORMAT(5HO J =912910Xs8HUN(J) =9F15e3s10Xs11HUNNEW(J) => 
петрова 
UN(J) = UNNEW(J) 
КЕЧКЕ ໆ ç H _— — و‎ í—————— — ИИ 
GO TO 101 


107 WRITE (6516) J»U(J) 
16 FORMAT (5HO J =912910X98H U(JU) =»F15.3) 
101 CONTINUE 
WRITE (69151) Nl» N2» МЗ 
151 FORMAT (/36X93A6/) 
WRITE (69901) 
901 FORMAT (/12X»8HP(J)DATA+912X»BHU(J)DATA»+11xX9+9HPX(J)DATA » 
111X»9HUN(J)DATA/) 
19HUN( J)DATA/) 
WRITE (691000) (P(I)s U(lI)s PX(1) > UN( I)» I = lM) 
1000 FORMAT (4(5X»F15.3)) = = E 
WRITE (653001) 
3001 FORMAT (/12X»BHP(J)DATA»>12X»9HPR(J)DATA>+10X>9HPX(J)DATA» 
110X»10HPRN(J)DATA/ ) 
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192 


1100 


WRITE (691000) (P(I)»>PR(I)»>PX(TI)>PRN(I)>9 I =1»M) 
FORMAT (F1005) 

FORMAT (12) 

WRITE (79152) Nls N2s N39 N6 
FORMAT (36X93A6920X 9A6 ) 

WRITE (791100) (UlK)» K =ls M) 
WRITE (791100)(UN(K)9 K 71» M) 
FORMAT (1X96(F10e493X) ) 

WRITE (692) 

FORMAT (1H1) 

6O TOLI 

END 
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FORTRAN PROGRAM FOR ORIGINAL LOGARITHM METHOD 
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Appendix V 


COMPUTER PROGRAM TO CALCULATE 
COEFFICIENTS FOR A UTILITY CURVE 


Purpose 


The purpose of this program is to calculate the coef- 
ficients for a utility curve from the utile-dollar data 


points computed by one of the data smoothing programs. 


Language 
Fortran IV (IBM 7040 Computer). 


Symbolic Dictionary 
VARI ABLE S/A* I/O** DESCRIPTION 


N o I Total number of data points. 

U A I&O Utility values for the data points. 

X A I&O Dollar values for the data points. 

B 5 0 The unknown variable within the log- 
arithm function which causes the 
curve to go through the (0,0) point. 

UCAL A O Calculated utility value for tne 
data points. Computed internally 
with the use of the computed coef- 
ficients. 

DEVIAT A O Deviation of data point utility and 
the calculated utility values. 

STDEV S O Standard deviation of the data point 
utility and the calculated utility 
values, 

A 5 Q The unknown variable which controls 


the magnitude of the curve. This is 
computed internally. 





*S - Single variable; A - Array of variables. 
**TI = Input; O = Output. 






BOB * ALOG (BOLD) 


FOLD E De EUN 54. , 108 ) “2 (ຖ້ , 108 ) ) 
-Q (N, SAM, U, BOI 
"ກ 


20 


3AM(J) = ALOO (X(J) + BOLD) 


rLOW DIAURAM FOR 


ONIQINAL LOJANI THA ИКООНАМ 















ОНПКИКО DOLLAH AND 
UTILITY VALUES CAN 
HE HNAD IN FHOM 

Е: ТНИН ОК ТИЕ DATA 
SMOOTHING PHOORAMS 


=» у 
CALL CLOCK (J) 


‚М = KMAX * MMAX * ) 
HOLD *-X(1)«91,001 
HTHY * -X(1)91.002 
ЕР = ¿005 

B = BOLD 

NCIEK * 1 

NCTH = O 


9e (X,N,B,SAM) = B(X,N,B,900) 





XJB = X(J) * B 


:: р 
= 


НС КЕК 






A 
NCHEK = NCIEK * 1 
B = (BOLD + B)/2 


00 TO 200 





X (K 
U(K 





* X(KX 
- || 











HOH = ALO (В) 


=ч (Н, ЗАМ, О ap rdi 
DELL = Porne (b - BOLD іе ~ РОРВ) 
В = в + DELbe1,1 + „OUS@DELB/ABS(DELB) 









— 


NCTH " NCTM * 1 


<s 


(ғ) OALCULATIONS 
DID NOT COMVEROR 





WRITE 
UTILITY CURVE PIT 


WRITE 
A, Al, B 


| — WRITE 
XDATA, UDATA, UCALC, DEVIATION 


UCAL(J) * Ae(SAM(J) - воз) 
DEVIAT(J) = U(J) < UCAL(J 







NHITE 
X(J3), U(J), UCAL(1), DEVIAT(J) 


SUM = 0,0 


SUM=SUM*(DEVIAT(I))eDEVIAT(I)) 


30 


AM = Н „ | 
STDEV=SQrT (SUM /AM ) 


WHITE 
STDEV 
J 
2 CALL CLOCK(J) 


<. +) 


РОТ (Хх, м, в, И) (# (Н, ЗАМ) -У(Н, ЗАМ, 300) 92 (ш,в 
Н, В, ЗАМ ) -3(Х,Ы, 8,808 


WHITE 


Ди 


4 
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Су Гу су CY €CY €Y CY CY CY CY €Y C боб 


THIS IS THE ORIGINAL LOGARITHM PROGRAM. 


PROGRAM DESIGNED TO FIT AN EQUATION TO UTILITY CURVES 
BASIC EQUATION. EPS - U(J) — A*ALOG(X(J) *B) - A*ALOG(B) 
TO SIMPLY THE FORMULAS», THE FOLLOWING TERMS ARE SET 
EQUAL TO REOCCURRING GROUPS OF VALUES, 
BOB < ALOG(B) 
SAM(J) = ALOG (X(J) + B) 
Н = ((1./X(J)+B)-10./B) 
DELEPS/DELAZA*SUM( SAM(J)*SAM(J) ) -A*X2, * SUM( SAM(J) *BOB) * 
AX*SUM(BOB*BOB ) - SUM ( SAM-BOB ) 
DELEPS/DELB =A*SUM((SAM(J)-BOB)*H) — SUM(U(J)*H) 
TO SIMPLY CALCULATIONS» THE FOLLOWING TERMS ARE SET 
EQUAL TO REOCCURRING GROUPS OF VALUES. 
ULCJ) *XCSAMUJ) -BOB) 
SAM ( J) *H 
BOB *H 
U(J) *H 
2e *SAM(J) *BOB 
SAM( J) ¥¥2 
BOB*BOB 
T/(R—S) = Q/(W-V+Z) 
DIMENSION X(20)» ЗАМ( 20)» (20) » ЦСА| (20) » DEVIAT (20) 
2 READ (5910) BTRYs BOLDs No EP 
READ (5320) (X(J) s U ())» J = JN) 
WRITE (692000) 
WRITE (692100) BTRY» BOLD» N» EP 
2000 FORMAT (7Xs4HBTRYs15X»4HBOLD»15Xs IHN»15X,» 2HEP///) 
2100 FORMAT 14X»>F10.3s9X»F100e3s10X»12»14X»F6.4///) 
WRITE (651111) 
1111 FORMAT (13X>1HB>15X»4HDELB/7/) 
10 FORMAT (2Ғ10.3»912» Ғб.й) 
20 FORMAT (2F10.3) — 
B=BOLD 
NCHEK = 1 
NCTR = 0 
DO 100 J= 198 
100 SAM(J) = ALOG (X(J) + BOLD) 
BOB = ALOG (BOLD) 


PNZ HAN DO 


FOLD = T(X»N»B»U) * ( W(N9SAM) — V(N>SAM9BOB) + Z{NsBOB) ) 
1 -Q(N>SAM+9U»+BOB)*R(X>N9IB9SAM) — S(X9N9B9BOB) 
B = BTRY 


180 DO 300 J = 1» М 
200 XJB = X(J) * B 
IF (XJB) 210, 210, 300 
210 60 ТО ( 220» 220» 220» 290) » МСНЕК 
220 NCHEK = NCHEK + 1 
B = (BOLD + B)/2e 


Ж 









7. ., MTP OWT A ເ... , 3 


EET е mide d: ics 


баба у '19 Я! 09491220 TAA e 
`В - AT = ah „мо ram на 
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КЛЕ = ім 
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и ін ctt oni ໄ 
"ກາມ . ретро 1645 hr m Í z 

р = ұғ + +. | 

1. hb e 34) в € 


| lnot p 








290 


300 


9981 


400 
450 
3000 
500 


600 


30 


31 
700 
800 
850 


900 


GOTO 200 

B = - X(1)*1.001 

GO TO 500 

SAM(J) = ALOG (XJB) 

BOB = ALOG (B) 

FOFB - T(X»sN»B»UD) * ( W(N»SAM) - V(N»SAMS, BOB) 


1 -QIN»SAM»U»BOB)I#R(X»N»B»SAM) - S(X9N9BsBOB) 


DELS =. FOF =“ ВО N ЦИ OB OHIO) 
В = В + DELB * lel + e005 * DELB/ABS(DELB) 
WRITE (6,9981) B » DELB 


ҒОКМАТ (7Х>Е10.3%910Х»Е10.3) 
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DEVIAT(J) = U(J) - UCAL(J) 
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SUM = SUM + (DEVIAT(J))*(DEVIAT(J)) 

AM = N - 1 

STDEV = SQRT (SUM/AM) 

WRITE (6,31) STDEV 

GO TO 2 

FORMAT (21HOSTANDARD DEVIATION =9F10-5) 
FORMAT (1Н1»31Х918Н UTILITY CURVE FIT///) 
FORMAT (7X>3H A=>F10+.3+10X>3H B=>F10.3//) 
FORMAT (11X»5HXDATA» 15X55HUDATA» 15X»S5HUCALC» 
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FUNCTION Q ( 

DIMENSION SAM 

Q = 0.0 

DO 10 JF 1 
Q = Q 

RETURN 

END 
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FUNCTION R ( X +> No Bs SAM) 
DIMENSION X (20) + SAM (20) 
R = 0.0 
DO 10 J = 1» N 
К = R + SAM(J) * ( 1. 
RETURN 
END 
FUNCTION T ( Xs N+ Bs U) 
DIMENSION X (20)» U (20) 
T = 0,0 
DO 10 Jj = 1» М 

Ts T PUJ) # TAX 
RETURN 
END 
FUNCTION W ( N» SAM) 
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APPENDIX VI 
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